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ABSTRACT 

Studies of deep photometry of galaxies have presented discoveries of excess light in surface-brightness and colour profiles at large 
radii in the form of diffuse faint haloes and thick discs. In a majority of the cases, it has seemed necessary to use exotic stellar 
populations or alternative physical solutions to explain the excess. Few studies have carefully scrutinized the role of scattered light 
in this context. 1 explore the influence of scattered light on ground-based observations of haloes and thick discs around edge-on 
galaxies, haloes around face-on disc galaxies, host galaxies around blue compact galaxies (BCGs), and haloes around elliptical galaxies. 
Surface-brightness structures of all considered types of galaxies are modelled and analysed to compare scattered-light haloes and 
thick discs with measurements. 1 simulate the influence of scattered light and accurate sky subtraction on simplified Sersic-type and 
face-on disc galaxy models. All galaxy models are convolved with both lower-limit and brighter point spread functions (PSFs); for a 
few galaxies it was possible to use dedicated PSFs. The results show bright scattered-light haloes and high amounts of red excess at 
large radii and faint surface brightnesses for nearly all types of galaxies; exceptions are the largest elliptical-type galaxies where the 
influence of scattered light is smaller. Studies have underestimated the role of scattered light to explain their surface-brightness profiles. 
My analysis shows surface-brightness profiles that include scattered light that are very similar to and overlap measurements at all radii. 
The derivation of physical properties of haloes, thick discs, and BCG hosts from diffuse dafa is misleading since accurate and radially 
extended PSFs are non-existent. Significantly improved analyses that include new measurements of PSFs are required to study diffuse 
haloes furfher. 
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1. Introduction 


Diffuse excess light has been discovered around objects as diverse 
as galaxies, supernova remnants, planetary nebulae, and stars. For 
galaxies, it has been found around, for example, edge-on and 
face-on disc galaxies, low surface-brightness galaxies (LSBGs), 
blue compact galaxies (BCGs) where it is thought to be part of 
the host component, and elliptical galaxies. Excess light seen atop 
a single-component intensity profile near the centre region of the 
galaxy is referred to as a thick disc or a halo at larger radii. Stellar 
populations of haloes are determined through studies of radial 
colour gradients, which often show strong discrepant red excess; 
the shape of the surface-brightness profile and the halo thereby 
provide valuable Information about the galaxy structure and evo¬ 
lution. Haloes are also seen around galaxies in observations of 
resolved stars, but such observations are by necessity limited to 
more nearby galaxies. Haloes are typically very faint, and there¬ 
fore difficult to observe; because of this, the discoveries are for 
some types of objects limited to a small number. However, the 
number of new discoveries of haloes can be expected to increase 
as the number of deep observations increase (for example, van 
Dokkum et al. 2014; D’Souza et al. 2014; Due et al. 2015; Zheng 
et al. 2014). 


* Figures E.4 and E.5 are only available in electronic form via 
http://www.edpsciences.org 

** E-mail: CSandin@aip.de 


Few haloes have been found around edge-on disc galaxies. 
Excess light has more often been observed closer to the cen¬ 
tre of the disc in so-called thick discs, which presence is even 
reported as ubiquitous (Yoachim & Dalcanton 2006, hereafter 
YD06). Surface-brightness structures around face-on disc galax¬ 
ies sometimes show flatter intensities at larger radii, beginning at 
a so-called anti-truncation radius; such profiles have even been 
assigned their own class, type Iff (Erwin et al. 2005; Pohlen & 
Trujillo 2006; Erwin et al. 2008, hereafter EPB08), extending the 
two previous types first defined by Freeman (1970). In BCGs, the 
occurrence of star formation in the centre regions is determined 
by properties of the surrounding stellar host. Structural properties 
of hosts of BCGs give clues to the gravitational potential where 
starbursts take place. Although, the host galaxy is typically much 
fainter than the starbursts and is therefore difficult to observe. 


Observations of haloes and thick discs are challenging and 
require careful assessments of data-analysis issues such as the 
accuracy of flat fields, sky subtraction, and other data-reduction 
properties. Accurate sky subtraction is essential with faint galax¬ 
ies where haloes are visible at intensities a fraction of the intensity 
of the night sky. The influence of scattered light from the instru¬ 
ment, the telescope, and the sky has not been checked as carefully 
as the sky subtraction. Most studies merely conclude that their 
intensity profiles are shallower than the point spread function 
(PSF), which defines and measures the extent of scattered light, 
and that their data therefore are unaffected by it. 
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In Sandin (2014, hereafter Paper I), I overview the existing 
sample of radially extended PSFs, and conclude that the PSF has 
rarely been measnred ont to large radii (angles) and is mostly 
inaccnrate. The paper describes how the PSF can be measnred 
more easily at large radii, and that it is necessary to acconnt for 
variations of the PSF with time and wavelength, as well as spa¬ 
tially away from the optical axis. Example models of edge-on 
disc galaxies illustrate that it is necessary to analyse observations 
of edge-on galaxies with PSFs at least 1.5 times as extended as 
the measurements, or the influence of scattered light is nnderes- 
timated. An alternative solntion regarding the origin of the halo 
around NGC 5907 is also presented. My re-assessment of the ac¬ 
curacies of the sky and the PSF suggests that the halo observations 
of NGC 5907 are mainly explained as scattered light. 

The mathematically correct procednre to remove scattered 
light is to deconvolve the measured image with the corresponding 
PSF. This is diflicnlt, however, as a deconvolntion reqnires a 
high accuracy at all intensities, both in the data and in the PSFs. 
Another approach is to probe the inflnence of scattered light by 
modelling snrface-brightness profiles that are at first convolved 
with the PSF and thereafter compared with the measnrements. 
Notably, such tests have been made for elliptical galaxies (for 
example, Capaccioli & de Vanconlenrs 1983; Michard 2002, 
hereafter M02) and edge-on disc galaxies (de Jong 2008, hereafter 
108). The studies on elliptical galaxies And that effects of scattered 
light are small. Some 20-80 per cent of the excess light in data 
of the Sloan Digital Sky Survey (SDSS) has, initially, been found 
to be scattered light (108). 

The qnestion is how scattered light more generally affects 
the faint onter parts of snrface-brightness profiles of galaxies 
observed in diffuse light. If scattered light is a key effect, it implies 
that conclusions of observations need to be reconsidered with a 
more carefnl analysis that involves accnrately measnred PSFs. 
My simnlation-based study scrutinizes haloes and thick discs 
in galaxies of the different kinds mentioned above to show that 
scattered light indeed is such a key effect. I largely use the same 
PSFs as in Paper I, which span lower-limit to larger effects for all 
kinds of instrnments and telescopes. The modelled profiles are 
compared with available measurements, and excess is studied in 
one colour for each galaxy. 

As in Paper I, this study focuses on ground-based observations 
in the visnal wavelength range, 300-900nm. The focns is also 
on symmetric structures, and I ignore mnch fainter stmctnres 
that show parts of shells or tidal tails and streams as are seen 
in, for example, Martmez-Delgado et al. (2010) and Dnc et al. 
(2015). The analysis method used to model all intensity profiles 
is described in Sect. 2. 

I scrntinize analyses of observations of different types of 
galaxies in separate sections. The case studies are ordered accord¬ 
ing to how strongly observations of the specific types are affected 
by scattered light, beginning with the most affected objects; most 
stndies have delimited themselves to one galaxy-structnre type 
following the same scheme. Three additional haloes around edge- 
on galaxies are studied in Sect. 3, which also contains a discussion 
of the thick discs aronnd edge-on disc galaxies of YD06, and a 
critical stndy of the finding of Bergvall ef al. (2010, hereafter 
BZCIO), who conclnde fhaf red excess aronnd edge-on LSBGs 
is not explained by scattered light. Five face-on disc galaxies are 
studied in Sect. 4, with a focus on the type III profile. Numerous 
surface-brightness profiles of galaxies described with a Sersic 
function are discussed in Sect. 5, including five BCGs and ellipti¬ 
cal galaxies. Each section begins with an overview of studies of 
the respective type of galaxy structure to indicate the extent and 
potential impact of this work. Several of the studied galaxies are 


discussed in appendices to corroborate the results in the main text. 
Section 6 contains a brief discussion including consequences for 
futnre observations, and the paper is closed with conclnsions in 
Sect. 7. 


2. Method and its application to example models 

The method of the simnlations and the analysis of each galaxy 
is split in three parts, as in Paper I. In Sect. 2.1, I select a set 
of PSEs nsed to estimate varying scattered-light effects. A set 
of models of the surface-brightness structnre is conflgnred in 
Sect. 2.2. Section 2.3 describes how the PSEs are applied to 
the model strnctures and the ontcome is analysed. In addition, I 
illustrate how the radial extent of the PSFs and the sky-subtraction 
accnracy affect models of the considered types of galaxies in 
Sect. 2.4. 

2.1. Choosing representative PSFs 

The PSFs used in this study should describe both temporal varia¬ 
tions and the red-halo effect in the i band (as described by Sirianni 
et al. 1998), and they should extend out to r ^ 900". As in Paper I, 
I chose to use the only published extended PSFs of M02, which 
were measured with both the (Cousins) V and the (Gunn) i bands, 
at two distinct occasions separated by three months: PSFvom and 
PSF, om were measured three months before PSFv 3 ni and PSF,- 3 n,. 
The difference with radins of PSFy^m-PSFyom on average is 
higher than that of PSF, 3 n,-PSF, om (see Fig. A.l). The radially 
extended PSFk 7 i of King (1971), which was measured (mostly) in 
the B band is used as a general lower-limit indicator of scattered 
light (cf. Paper I). The comparison of the PSFs in figs. 1 and 2 in 
Paper I demonstrate that PSFyom is ^ 0.5 mag arcsec^^ brighter 
at intermediate radii (20 ^ r ^ 60") than the other treated PSFs. 
PSFv 3 ni, meanwhile, is another 0.5 magarcsec^^ brighter. In Pa¬ 
per I, I reassessed PSFmbh that was measured by Morrison et al. 
(1994) in their study of NGC 5907; the new PSF“"jj lies between 
PSFyom and PSFyom ■ PSFj Om appears to represent the average 
/-band SDSS PSFs well where r k, 10". Some PSFs show more 
light at large radii than the r^^ power-law of PSFkti; the predic¬ 
tive ability of models of larger objects is weaker, as the far wings 
of the PSF are generally poorly determined. The complementary 
PSFs discnssed in this paper, which all extend beyond r ^ 10", 
are shown in Fig. A.l. 

In the simnlations, I assnmed that the two V-band PSFs, as 
well as PSFk 7 i, are the same in the B, g, R, and r bands, as well as 
the photographic J band (that is measured using a Illa-J emulsion 
and a Wratten-2C Alter), and that the two /-band PSFs are the 
same in the I band; the colour predictability is hereby delimited 
to one colour, for example, V - i. I do not differentiate between 
Consins, Johnson, or other photometric systems. All five PSFs 
were measnred af a seeing of several arc seconds, which is why 
the spatial resolution in the brightest centre region is poor. Here, 
faint diffnse emission is stndied, where the spatial resolntion 
and resulting lower intensities in the centre regions are of minor 
importance. 

2.2. Configuring sets of modet surface-brightness structures 

I discriminate between galaxies, which surface-brightness struc¬ 
tures are fitted well with models that only vary with the radins, 
and disc galaxies viewed at low (face-on) and high (edge-on) in¬ 
clination angles. The relation ju - -2.5 log[o(/) is used to convert 
between magnitudes p and intensities I. 
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The common Sersic function (Sersic 1968) is used to model 
surface-brightness structures of, for example, elliptical and lentic¬ 
ular galaxies, cD galaxies, bulges of disc galaxies, and BCGs, 


where r is the radius, n the Sersic-profile index, rg an effective 
radius that encloses half of the total light of the profile, /g an ef¬ 
fective surface brightness at rg, and bn -2n- 0.327 (for example, 
Capaccioli 1989; Caon et al. 1993). Also, 7(0) = /g exp (bn)- For 
exponential profiles (n - 1 .0), the function can be replaced with 
a representation that instead uses a scale length h, 

/(r) =/oexp^-^j, (2) 

where Iq - /gexp(7>i) is the surface brightness at the centre, 
and h - r^fbi ^ 0.598rg. Moreover, when modelling flattened 
spheroidals, it is assumed that the minor-to-maj or axis intensity 
ratio equals the c/a ratio. The ellipticity (flattening) is defined as 
e = 1 - c/a. Otherwise, any azimuthal dependence of measured 
profiles is ignored. Sersic-profile surface brightnesses reach very 
high values at sub-arc second angles for higher values on n. To 
scale the models reasonably, each model Sersic profile is normal¬ 
ized with the model centre intensity that results after convolving 
the model with a 1" seeing profile; I make a note of resulting 
non-zero offsets in figures of affected profiles. 

Disc galaxies are suitably described in cylindrical coordinates; 
assuming an isothermal disc, the space-luminosity density can be 
described by (van der Kruit 1988; van der Kruit & Searle 1981, 
hereafter KS81a) 

I'(r’,z') = /o oexp^-^j x 2^2^"^sech2^"* , (3) 

where 7^ ^ is the centre intensity, r' the radius, the scale length, 
z' the vertical distance from the centre, zo the vertical scale height, 
and Hs is set to 1, 2, or oo. The intensity drops to zero at a galaxy- 
specific truncation radius. When a disc galaxy is projected edge 
on, there is a(n optional) Sersic-profile bulge, and the truncation 
radius and dust extinction are ignored, the surface-brightness 
structure becomes 

I(r,z) - 7(r) + 70.0^77, x 2-2/«^sech2/”* , (4) 

where 7o_o = 2/1^7^ ^ is the centre intensity, r the major-axis ra¬ 
dius, z the vertical-axis distance from the major axis, and the 
modified Bessel function of fhe second kind. The intensify is 
slighfly lower al larger radii when fhe truncation radius is finite 
and taken into account; considering how uncertain the PSFs are, 
this effect is ignored here. The vertical structure can be assumed 
as ‘isothermal’ (n^ — 1), exponential (wj = o°), or in between 
(Hs = 2). When Wj = oo, the z-dependent function reduces to 
exp(-z/zo); papers that discuss exponential structures often use 
the notation — zo- No attempt is made to make a perfect fit of 
the centre regions of edge-on galaxies. 

Finally, with disc galaxies that are instead viewed face on, 
a surface-brightness structure is used that accounts for radial 
breaks in the intensity structure as well as a Sersic-profile bulge. 
Following the approach of Bakos & Trujillo (2012), 

7(r, z) = 7(r) -H 7o,i exp ^ j ® ('') 

-H 7o,oexp^-^j(l-©(r)), (5) 


where Iq - IzqIq q is the centre intensity, 7o,i and Iq o (hi and ho) 
are the centre intensities (scale lengths) of the inner and the outer 
discs, respectively, and the step function 0 (r) = 1 if r < rbrk (rbrk 
is the break radius), and otherwise zero. 

2.3. Models and measurements analysis procedure 

Each two-dimensional V-band (as well as B-, g-, R-, r-, and 
7-band) model image is convolved individually with the three 
resampled and normalized two-dimensional V-band PSFi/om and 
PSFj/ 3 n,, and the B-band PSFkvi. This is repeated for each i- 
band (as well as 7-band) model image with PSFkOm and PSF; 3 m- 
Unless specified otherwise, the /-band models use the same model 
parameters as the V-band model, with two exceptions; r^j - 
0.95re_y (hi - 0.95hv) andzo.i = 0.95zo,v, which result in negative 
slopes in y - / with increasing radius. Thereby, any red-excess 
haloes in convolved models are induced by the PSF. Except the 
LSBG models in Sect. 3.3,1 made no attempt to fit observed V-i 
colours, but po,v-Pe,i, Po,v-PQ,i, and /Uo,o,v-po,o,i are set to values 
that result in a very rough general agreement between models and 
measurements; for example, I used pv - Pi = 1.3 mag arcsec^^ 
with all edge-on disc galaxy models. The used PSE image is 
always twice as large as the model image, to avoid PSF truncation 
effects in convolved images (see the outcome of the toy models for 
edge-on galaxies in Paper I and for Sersic-type and face-on disc 
galaxies here, in Appendix B). The PSF and the model images are 
resampled to use the same pitch and about 100-200 pixels on the 
side, typically, to keep calculation times short. All convolutions 
are made by direct integration. 

Convolved surface-brightness profiles are plotted together 
with model surface-brightness profiles for a cut in the image that 
extends from the centre and outwards. For edge-on galaxies, this 
cut is always positioned on the minor axis, directed outwards 
from the centre on the same (vertical) axis. In cases where deep 
enough B-, V-, g-, R-, r-, or 7-band measurements exist in the 
literature, they are plotted as well; where tables were missing, 
I used dexter' (Demleitner et al. 2001) to extract the data. The 
three resulting V — i colour profiles are shown in a separate lower 
panel; the model, convolved models that use the earlier PSFyom 
and PSF,_om, or the later PSFy^m and PSF, 3 m- 

The PSF induces scattered light, where surface-brightness 
profiles of input models and convolved models differ. Larger dif¬ 
ferences between profiles of convolved models that use PSFyom, 
PSFv 3 m, and PSFkti in the V-band, or PSF,;om and PSFj^m in the 
/-band, are indicative of stronger time dependence in the PSFs. 
The opposite applies when such differences are smaller. Mea¬ 
sured values due to scattered light should fall atop profiles of a 
model convolved wifh fhe PSF al fhe time of fhe observations. 
For some telescope and instrument setups, the PSF contains even 
more energy at large radii than PSFv 3 m, and in this case it is 
necessary to convolve affected models using an even brighter 
PSF. The PSF is mostly not measured or published, and in this 
case measurements expected to be scattered light plausibly fall 
between profiles of convolved models that assume the lowest and 
highest amount of such scattered light (PSFkvi and PSFy^m). The 
predictive power is weaker without a published PSF, but models 
nevertheless provide valuable hints of the PSF dependence of 
specific model parameter setups. Finally, measurements that are 
brighter than the convolved model that uses the correct PSF are 
real and are not (only) due to scattered light, and measured values 
that are fainter than the same convolved model are non-physical. 

' DEXTER is available at: http://dexter.sourceforge.net. 
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Table 1. Results of the models of the sky, using the intermediate-size 
configurations of table 3 in Paper I and Table A.2. Values in rows printed 
in boldface mark input models and remaining values models convolved 
with PSFv.om- 


Model 

ds 

d = 

= 0.05 

d = 

= 0.10 

d = 

= 0.15 

Him 

Asky 

Him 

Asky 

Him 

Asky 

edge-on 

7.5 

18.0 

4.4^5 

19.0 

4.9-5.1 

20.0 

5.4-5.8 


7.5 

36.0 

4.3 

- 

- 

- 

- 


5.0 

15.0 

1.8-1.9 

18.0 

2.5-2.8 

20.0 

2.9-3.2 

Sersic 

7.5 

18.2 

4.2^4 

19.3 

4.9-5.1 

20.4 

5.1-5.5 


7.5 

20.0 

4.2^.3 

22.0 

4.9^.9 

24.0 

5.2-5.4 


5.0 

13.0 

1.6-1.7 

15.0 

2.4-2.6 

16.0 

2.7-3.0 

face-on 

7.5 

22.9 

3.5-3.6 

24.2 

4.3-47 

25.7 

4.7-5.0 


7.5 

27.0 

3.5-3.6 

32.0 

4.2^.3 

38.0 

4.6^.8 


5.0 

18.0 

1.0-1.3 

20.0 

1.8-1.9 

22.0 

2.2-2.6 


Notes. Column 1, intermediate-size galaxy model; Col. 2, error of added 
and subtracted sky (mag arcsec^^); Cols. 3 and 4, radius (arcsec) and 
magnitude range (mag arcsec^^) at the innermost radius where the differ¬ 
ence between the model and the model with added and subtracted sky 
is d ^ 0.05 mag arcsec^^; Cols. 5 and 6 (Cols. 7 and 8), same as Cols. 3 
and4 ford - 0.10magarcsec^^ (d - 0.15magarcsec^^). 

A scattered-light halo radius ri m is defined as the lower limit¬ 
ing radius where the convolved model intensity at all larger radii 
is > 10 per cent higher than the model intensity. Note that rno 
depends on both object parameters and the PSF. The usefulness 
of riio at larger radii (r ^ 100") is at the moment doubtful in 
comparisons with real galaxies, since PSFs are nearly unknown 
there. 

2.4. Using example models to test the influence of the radial 
extent of PSFs and the accuracy of the sky subtraction 

I demonstrate effects of the use of truncated PSFs in the analysis 
of observations of edge-on galaxies in Paper I. In Appendix B, I 
extend the study with completely analogous models for face-on 
disc galaxies and galaxy structures described with Sersic-type 
profiles (where n = 1). The results of the demonstration are 
unambiguous - surface-brightness structures of these galaxies 
are also strongly affected by scattered light, and effects are only 
somewhat dependent on the size of the galaxy. These new models 
indicate that the required extent of the PSFs for these structures is 
some 1.1 times the radius of the measurements. Depending on the 
brightness of the PSF, the halo is up to 2.5 mag arcsec^^ brighter 
around the edge-on galaxy models than around the Sersic-type 
galaxy models. The face-on galaxies are an intermediate case 
where the brightness of the halo depends on the properties and 
extent of the galaxy inside a (possible) break radius. I show in 
Sect. 5 that effects are smaller in elliptical galaxies (n = 4) and 
cD galaxies (n ^ 8). 

Accurate subtraction of the sky is imperative, considering 
that faint haloes and thick discs are observed around already 
faint galaxies. It is beyond the scope of this paper to provide a 
general study of the influence of sky subtraction for all possible 
configurations of galaxy models and levels of the sky background. 
Instead, I show three examples of the already modelled galaxies 
to illustrate the interplay between sky subtraction and scattered 
light. The three models are the intermediate-size models of an 
edge-on disc (Paper I), a Sersic-type, and a face-on disc galaxy. 

I simulated the influence of the sky by adding a constant 
value to each modelled two-dimensional structure (arbitrarily) 


3.8magarcsec“^ fainter than the centre intensity. Thereafter, I 
subtracted two constant values in separate models offset from 
the added value by +1 per cent (ds = 5.0magarcsec“^) and 
+0.1 per cent (Ss -7.5 mag arcsec^^), respectively. The added- 
and-subtracted sky models were then convolved with PSFv;om- 
Resulting profiles are shown for all three galaxy types in Fig. 1. 
I measured minimum radii and corresponding magnitudes 
where the sky-subtracted models deviate from the input model by 
6 > 0.05, 0.10, and 0.15 mag arcsec^^; these values are presented 
in Table 1. 

The plots clearly show that the accuracy of the sky-affected 
input models are several magnitudes lower than the added per¬ 
turbation of 6s - 7.5 mag arcsec^^ when the allowed deviation 
6 < 0.15 mag arcsec^^. The limiting magnitudes for the edge- 
on (Sersic-type) galaxy are then 4.4^.5 (5.4-5.8) mag arcsec^^ 
fainter than the sky when 6 = 0.05, 0.15 magarcsec^^; the cor¬ 
responding radii are r\i„i - 18, 20". These values correspond 
to a lowering of the accuracy of the sky by 3.0-3.1 (3.1-3.3) 
magarcsec^^ when 6 - 0.05 mag arcsec^^ and by 1.7-2.1 (2.0- 
2.4) magarcsec^^ when 6 = 0.15 mag arcsec^^. The values are 
slightly lower for the more extended face-on disc galaxy model, 
where the limiting magnitude is 3.5-3.6 (4.7-5.0) mag arcsec^^ 
for 6 - 0.05 (0.15) magarcsec^^, at the radius = 23" 
(riim = 26"); in this case, the accuracy of the sky is lowered by 
3.9^.0 (2.5-2.8) magarcsec^^. Surface-brightnesses measured 
at larger radii can be compared with models to show if they may 
be affected by scattered light. However, it is only possible to 
recover an exponential profile, as in this case, out to a maximum 
radius determined by the accuracy of the sky. 

Limiting magnitudes are only slightly lower in the convolved 
models than in the input models; the tabulated values show the 
range 0.0-0.4 mag arcsec^^, considering all three values of 6 for 
all three types of galaxies. The convolved sky-affected models 
that used a perturbation of ds = 5.0 mag arcsec^^ deviate from 
the convolved models at smaller radii and much brighter inten¬ 
sities; the values are very close to 2.5 mag arcsec^^ lower than 
the models that use a ds -7.5 mag arcsec^^ perturbation. For the 
edge-on galaxy model, and compared to the model that used a 
7.5 mag arcsec^^ perturbation, the magnitude is 2.5 mag arcsec^^ 
lower with d = 0.05, whilst the limiting radius is rum = 15". 
The corresponding values of the Sersic-type (face-on disc) galaxy 
model are 2.5-2.6 (2.4-2.5) magarcsec^^ lower and rum = 13- 
16" (18-22") for all three values on d. 

3. About analyses of observations of faint 
structures around edge-on galaxies 

Thick discs and haloes are observed in resolved star counts of 
nearby galaxies and in photometry of both near and far galaxies; 
here, I only consider diffuse photometric measurements. Resolved 
star counts are always real, regardless of if stars are attributed to 
a thick disc or an extended thin disc (Bovy et al. 2012) - discrete 
counts are unaffected by scattered-light effects. 

Tsikoudi (1977, 1979) presents the first detection of diffuse 
excess light around three bright and nearby lenticular SO galaxies. 
She uses observations with photographic plates in the B and V 
bands. Burstein (1979), who also uses photographic plates in 
the B band, coins the expression ‘thick disc’ in his study of five 
SO galaxies that all show excess light. Following studies switch 
from photographic plates to CCD measurements and report on 
both positive and negative identifications of diffuse thick discs 
in a number of galaxies and different bands. The selection of 
objects for this study is made in Sect. 3.1, which also outlines the 
remaining parts of this section. 
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z (arcmin) r (arcmin) r (arcmin) 

0.2 0.4 0.6 0.2 0.4 0.6 0.2 0.4 0.6 



z (arcsec) r (arcsec) r (arcsec) 

Fig. 1. Surface-brightness profiles of the three intermediate-size example models, which illustrate effects of the accuracy of the sky subtraction 
using PSFv-om- The three panels show the: a) edge-on disc galaxy (intensities are drawn versus the vertical distance z instead of the radius r), b) 
Sersic-type galaxy, and c) face-on disc galaxy. In each panel, the galaxy model profile is drawn with a thick black solid line, and the convolved 
model profile is drawn with a thin blue line. The horizontal dashed red lines indicate the used background level of the sky, and the purple and blue 
short horizontal lines on the right-hand side of each panel 0.01 and 0.001 times the sky level. The thick lines overplotted with white dotted lines 
indicate the resulting profiles when the sky is under or over subtracted by 0.1 per cent (black and blue lines) and 1 per cent (purple lines). 


3.1. On scattered light and details of the object selection 

Few of the numerous studies of detections in diffuse light test the 
influence of scattered light on their observations. In general, stud¬ 
ies of ground-based observations merely note that their objects 
are unaffected by scattered light. Neither Elmegreen (& Elmegreen 
(2006) nor Comeron et al. (2011) address the effects of the ex¬ 
tended PSE in the analyses of their space-based observations. 108 
reports that edge-on galaxy observations are affected by scattered 
light, but does not specifically mention thick discs. In view of the 
results of Paper I, it seems that 108 underestimates the influence 
of the scattered light as he both neglects temporal variations of 
the PSE at large radii and extrapolates it with a too steep slope. 

Here, I extend the example of NGC 5907 in Paper I with 
three corroborating examples of diffuse haloes. An extended PSE 
was measured for the observations of the comparatively distant 
galaxy IC 5249, but its influence has not been checked. The thick 
disc and halo of NGC 4565 was measured multiple times; the 
resulting surface-brightness structures of the halo are suspiciously 
similar. No halo is found in analyses of diffuse observations of 
NGC 4244, but resolved-star counts do show a halo. These three 
galaxy haloes are scrutinized in the context of scattered light in 
Appendix C. To demonstrate that scattered light can give rise to 
diffuse fhick discs in small-to-intermediate size galaxies, I show 
two examples with the size-limiting galaxies FGC 310 and FGC 
1285 in Sect. 3.2. Finally, as a general attribute, the influence of 
the PSE has been dismissed in a study of stacked SDSS images of 
edge-on LSBGs (BZCIO). I show in Sect. 3.3 how this dismissal 
can be explained when accounting for temporally varying PSFs, 
and therefore scattered light cannot be excluded as an explanation 
to the surface-brightness structures for these galaxies either. De¬ 
tails of the observations and existing and new model parameters 
of all six considered galaxies are collected in Tables A.l and C.l, 
respectively. 

3.2. On the origin of the observed thick discs around the 
galaxies FGC 310 and FGC 1285 

Dalcanton Sr. Bernstein (2000, hereafter DBOO) and Dalcanton Sr 
Bernstein (2002) select a sample of 47 edge-on galaxies from the 
Flat Galaxy Catalog (FGC; Karachentsev et al. 1993). The authors 


carefully consider issues of calibration uncertainties, the influence 
of field stars, and sky subtraction. They measured PSFs for each 
of the observed bandpasses in the radial range 0 < r :$ 6", which 
they compare with the observed profiles, concluding that the PSE 
is too narrow to have any effect. 

These observations are worth a closer examination, since 
the authors omitted a calculation of the integrated light of the 
galaxies. DBOO present contour diagrams for each galaxy, where 
contours are drawn at a 1 mag resolution. Later, Dalcanton & 
Bernstein (2002) present minor-axis surface-brightness profiles, 
which, however, are not as easily interpreted (they appear to differ 
from the data of DBOO and YD06). YD06 present two model 
fits for each galaxy: one with a single disc and one where a 
thick disc is superposed on a thin disc. It is worth noting that the 
authors weight their two-dimensional fits to the pixels with the 
lowest signal-to-noise ratio; this approach could perhaps explain 
differences between the model profile fits and the profiles fhaf 
can be extracted from the contour diagrams of DBOO. I chose one 
of their smallest and the largest appearing galaxies FGC 310 and 
FGC 1285, respectively. The centre parts of the smaller galaxy 
are poorly modelled here since the seeing values of the PSFs I 
used are large (about 478, which results in somewhat lower centre 
intensities). 

I modelled FGC 310 with one disc where /i, = 874, zo - 176, 
and Hs = 1, Fig. 2a. The input model convolved with PSFkvi 
agrees with the single-disc model, and when it was instead con¬ 
volved with PSFvom h agrees with the superposed thin-and-thick- 
disc (two-disc) model. By chance, the input model agrees with 
the thin disc of YD06. The minor-axis profile that I extracted 
from the contour diagram in DBOO is more similar to the two-disc 
model, at least where r ^ 5". The colour profiles in Fig. 2b show 
red excess for z ^ 4", where A{R - i) ^ 0.44-0.50 mag arcsec^^ 
at z = 7" and about 0.92-1.2 mag arcsec^^ at z = 10". The lower 
(higher) values again used the latter PSEy^m and PSE/^m (earlier 
PSFyom and PSF, om)- Notably, the effects of scattered light are 
larger for the smallest appearing galaxy in the study of YD06, 
FGC 1063, because the model parameters are even lower than in 
FGC 310. 

I modelled the larger galaxy FGC 1285 with a bulge that is su¬ 
perposed on a disc. The parameters of the bulge are n - 1.0, 
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Fig. 2. Minor-axis surface-brightness profiles 
versus the vertical distance z for two edge-on 
disc galaxies, a) and b) FGC 310, and c) and 
d) FGC 1285. Blue and purple lines in panels 
a) and c) show f?-band profiles, and red lines 
t-band profiles. Model profiles are drawn with 
thick solid lines. Solid (dash-dotted) lines are 
profiles of convolved models using PSFv/om and 
PSF,_om (PSFv.sm and PSF, 3 n,), the purple line 
used PSFk 7 i. The black (grey) line shows the 
respective fitted combined thin-and-thick-disc 
model (single-disc model; of YD06). Crosses 
indicate f?-band measurements that were ex¬ 
tracted from fig. 3 in DBOO. The lower limiting 
radius riio - where the convolved models us¬ 
ing PSFi/oni, PSF,-,oni, and PSFk7i lie > 10 per 
cent above the input model - is marked with 
a coloured bullet that is surrounded by a black 
ring, b) Three colour profiles R - i are shown 
for: the model (thick solid line), the convolved 
model using PSTv-om and PSF, om (solid line), 
and the convolved model using PSFi/ 3 ni and 
PSF,,3„ (dash-dotted line). 


Te = 4"4, and — 19.3 mag arcsec^^, whilst the disc pa¬ 
rameters are = 17"3, zq - 9'.'6, po.o = 19.5 mag arcsec^^, 
and Hfi - 1. Minor-axis surface-brightness profiles are shown 
in Fig. 2c. For smaller radii r ^ 20", there are small differ¬ 
ences between profiles of the input models and the convolved 
models. The observations of this object are less strongly af¬ 
fected by scattered light than FGC 310. The single-disc pro¬ 
file of YD06 agrees poorly with the two-disc profile at all radii. 
In this case, when using PSFv/om, the scattered-light halo only 
becomes visible at an intensity that is some 4.8magarcsec“^ 
fainter than the centre intensity; with PSFkvi this value is about 

7.8 mag arcsec^^. The colour profiles in Fig. 2d show red excess 
forz ^ 15", where A(R- i) ^ 0.22-0.28 mag arcsec^^ atz = 23" 
and A(R- i) ^ 0.53-0.69 mag arcsec^^ at z = 30". 

The small sizes of all galaxies of YD06 imply that surface- 
brightness profiles and colours are highly sensitive to the PSF 
shape - and thereby to temporal variations of the PSF. The exact 
values of the model parameters that are fitted here are unimpor¬ 
tant. It should be clear that it appears to be of little value to fit 
a thin disc and a thick disc unless the data are first corrected for 
scattered light. Currently, the fitted components do not reveal any 
information about the real structure in the faint parts. To recover 
the correct structure out to z = 10", it is necessary to deconvolve 
observed structures with PSFs that extend to at least r - 15" in 
the considered bandpass (1.5 x 10", cf. Paper I). For FGC 310 
and using PSFkti (PSFyom), the convolved structure at z = 8" is 

1.9 (3.3) magarcsec^^ brighter than the model structure, which 
corresponds to a required S/N= 5.8 (21). The corresponding val¬ 
ues at z = 11" are 4.4 (5.8) mag arcsec^^ and S/N= 58 (210), and 
at z = 14" the values are 7.6 (9.0) mag arcsec^^ and S/N= 1100 
(4000), and the values increase for larger radii. The requirements 
are lower for FGC 1285, where the halo, however, appears at 
fainter intensities. At z = 40" the convolved structure is 0.56 
(1.5) mag arcsec^^ brighter than the model structure, the required 
S/N= 1.7 (4.0). 


3.3. A closer look at the role of time-dependent scattered light 
to the analysis of stacked Images 

Galaxy halo intensities typically reach levels much below one per 
cent of the background sky. One can assume that haloes share 
fundamental properties the same way galaxies do, to overcome 
the difficulty in measuring such faint intensities; van Dokkum 
et al. (2010) show that parameters of stacked images created from 
hundreds of synthetic galaxy models with random Sersic profiles 
mafch average values of sfacked profiles well. Fligh signal-to- 
noise values in halo measurements can then be achieved by stack¬ 
ing a large number of images of fairly similar objects. Zibetti 
et al. (2004), who are the first to use this approach, address the 
role of the PSF in stacked images of 1047 nearly edge-on disc 
galaxies of the SDSS; they find that haloes are an almost ubiqui¬ 
tous phenomenon around disc galaxies. The halo of their stacked 
image shows abnormal red excess that they attribute a physical 
origin, and not scattered light. J08 finds that the reason that the 
PSF plays no role is that Zibetti et al. (2004) do not consider 
the extended wings of the PSF properly - at least parts of the 
red excess in the halo can be explained by such extended wings. 
Despite accounting for extended wings of PSFs, 108 concludes 
that red haloes are still found far away from the disc (in both HST 
and SDSS images). 

BZCIO carefully analyse a stack of nearly edge-on LSBGs, 
and conclude that for this kind of galaxy the PSF is unimportant. 
I here demonstrate, based on the study of BZCIO, how both 
colours and faint stellar haloes can be induced as a consequence 
of temporally varying PSFs. 

BZCIO search for nearly edge-on LSBGs in g-, r-, and i- 
band SDSS images. They use the following four criteria to select 
galaxies: the mean surface-brightness inside the radius encom¬ 
passing half of the light in the r band must satisfy 23 < pr.so < 
25 magarcsec^^, the /-band isophotal diameter a > 10", the /- 
band axial ratio bja < 0.25, and the redshift z < 0.2. The halo 
region of interest that they discuss extends from 60 pixels and out- 
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wards on the vertical axis, which corresponds to radii r > 23'.'7. 
Their final sample contains 1510 galaxies that they analyse in 
four separate subsamples (A, B, C, and T = A + B + C), which are 
setup based on the g-r colour. Results are presented as statistical 
properties of each sample, where exponential scale lengths are 
measured in the range 1 - 16 kpc . The authors do not provide 
any information abont the identity or properties of individnal 
objects. They obtain one PSF for each of the three bandpasses 
from stacks of bright stars. 

Furthermore, BZCIO use two approaches to correct their re¬ 
duced images for scattered-light effects. In one approach, they 
deconvolve each image individually, before they stack them. In 
a second approach, they first fit a flattened-spheroidal and an 
exponential-disc model to a template stack of the 35 apparently 
largest galaxies. Thereafter, they scale the model to the same 
size as each galaxy in the respective sample and stack the scaled 
models. They provide fitted parameters for each bandpass of their 
template stack, as well as isophotal contonr plots of the stacked 
images of each bandpass and subsample. In their analysis of the 
stacked image data, they treat each snbsample individually in 
terms of a thin disc, a thick disc, and a halo. Their conclnsions 
are based on g - r and r - i colour profiles along fhe major 
axis for fhe fhin disc, and along fhe minor axis using wedges 
with an opening angle of 120° (thick disc) and 60° (halo). They 
find fhaf fhe r - i colour of fhe halo T sample decreases by 
r - i ^ 0.2 mag arcsec^^ due fo deconvolution. The remaining 
colour difference between more central regions and outer regions 
is A (r - i) ss 0.7 mag arcsec^^. Despite an individual deconvolu¬ 
tion of each image, they find red excess in the halo, and conclude 
that scattered light is not the reason for this excess. They instead 
explore a physical origin to explain the red excess. 

The authors do not consider temporally varying PSFs; the 
SDSS images are sampled during long periods of time, where 
such variations can be signihcant (cf. Paper I). Additionally, ob¬ 
jects of different geometrical properties are differently affected by 
the PSF - it is impossible to remove PSF effects exactly, without 
detailed knowledge of individual object properties and measure¬ 
ments. The exact PSFs are nnknown, and no attempt was made 
to reproduce the observations precisely. I used the following pro- 
cednre to simnlate the effects of ignoring temporally varying 
PSFs on the stack of galaxies. I convolved input g- and /-band 
surface-brightness models with the early-time PSFs (PSFyom and 
PSF, om), to give observed snrface-brightness structures. There¬ 
after, I deconvolved the result with the later-time PSFs (PSFy^m 
and PSF, 3 nj), to give PSF-affected input model structures. I used 
the Richardson-Lucy algorithm in the form of the freely avail¬ 
able program sgp for the Interactive Data Langnage (idl; Prato 
et al. 2012)^. I did not add any noise to the models or the PSFs, 
and using the same PSF I checked that the deconvolution always 
reproduced the input model structure when I deconvolved any 
convolved input model. I compared surface-brightness stmctures 
and the g-i colour of input models and PSF-affected input models 
with corresponding measurements. 

I hrst dehned the scattered-light affected input model. I mea¬ 
sured the isophote positions on the major and the minor axes 
of the g- and /-band contour plots of the B subsample (fig. 7 in 
BZCIO); the isophotes are separated by one magnitude. (I chose 
the B subsample as it appears less noisy than the other snbsam- 
ples in the same hgure.) Effects are expected to be slightly larger 
(smaller) for objects of the smaller (larger) objects of the A (C) 
subsample. I htted the major-axis values with two separate Sersic 
profiles where Ug - 1.0, r^ g — 26", p^^g — 24.0 mag arcsec^^, 

^ SGP is available at: http://www.unife.it/prisma/software. 


n, = 1.2, r^ i - 21", and pe,; = 22.3 mag arcsec^^. I then cre¬ 
ated two-dimensional surface-brightness stmctures by combining 
the two Sersic profiles with the major-to-minor axis flattened 
spheroidal parameters (cla)g - 0.34 and (c/a), = 0.50, which 
the anthors derive from the faintest isophotes in their template 
model stack. To have the PSF-affected input model reproduce 
the measurements on the minor axis in the isophotal image, it 
is necessary to nse lower valnes since scattered light affects the 
ellipticity at large radii (see the discussion for the BCG ESO 400- 
G043 in Sect. 5.2.1). I fonnd, by trial and error, that Ugj - 1.0, 
(cla)g - 0.20, and (c/a),- = 0.22 produce good results. Minor-axis 
plots of all models are shown in Eig. 3. 

Figure 3a shows a reasonable agreement between the two 
PSF-affected input models and the measnrements in both the g 
and the / bands; deviations are ;$ 0.25 mag arcsec^^. A scattered- 
light halo appears in the observed structnre at r ^ 12 " for the 
g band, and at r ^ 3" for the / band. In the PSF-affected input 
model, this radius is r ^ 12" for the / band. These scattered-light 
haloes appear with all four PSFs that I use, and they demonstrate 
that galaxy surface-brightness structures that can be represented 
with parameters such as these are strongly affected by scattered 
light already at small radii. 

The PSF-affected input model shows that the /-band scattered- 
light halo is more completely removed in the deconvolution 
than that of the g band, since at most radii (PSF, 3 ,^- PSF,_ 0 m) < 
(PSFy, 3 n,-PSFv,om)- Figure 3b shows that the colour prohle g - i 
of the input model decreases with the radius. A deviation from 
the model colour gradient - and a red-excess halo - results with 
the PSF-affected input model, where r S 16". The agreement 
between the PSF-affected input model and the g - i colour prohle 
of the deconvolved halo T sample is good where r ^ 15" (this 
was achieved by co-adding the two cnrves in flg. 14 of BZCIO, 
no similar prohle is presented for the B snbsample); deviations 
are $ 0.12 mag arcsec^^. The offset of r - / ^ 0.3 mag arcsec^^ 
where r 12 " conld arise due to the less excellent model ht in 
the g band (see the value at r ^ 7"), differences befween fhe B 
and the T samples, or due to PSF differences in the inner region, 
or because of all these reasons. 

My simulations cannot prove that the stacked-images halo 
arises only dne to scattered light, or that there is no red excess 
in the halo. However, the results provide strong indications of 
overlooked systematic properties that alone can explain both these 
features. Additionally, Tal & van Dokknm (2011) stack 42 000 
SDSS images of luminous red galaxies, and also use temporally 
averaged PSFs. These authors hnd support to the conclusions of 
Zibetti et al. (2004), in terms of decreasing r - / colour prohles in 
the inner parts of the stacked galaxy. The geometrical properties 
of their stacked image resnlt in models that are less affected by 
PSF effects than the stack studied here. Although, the outer parts 
of a model that I calculated for such luminous red galaxies show 
signihcant red excess also in this case (colours are not shown for 
this region in their paper). My analysis suggests that the complete 
neglect of extended PSFs in stacked images (for example, Zibetti 
et al. 2004, 2005; D’Souza et al. 2014) affects results more than 
just neglecting the temporal variations of the extended PSFs. 


4. About analyses of observations of faint 
structures around face-on disc galaxies 

Surface-brightness prohles of face-on disc galaxies are catego¬ 
rized into three types. Type I prohles are well described by a 
single exponential, whilst type II prohles show a break outside 
the main region of spiral-arm activity, which is followed by a 
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Fig. 3. Minor-axis g- and i-band surface- 
brightness profiles versus the vertical distance 
z of the LSBG model. Panels and lines are as 
described in Fig. 2. Circles (crosses) indicate 
measured values of the g band {i band; BZCIO). 
In panel a, the blue and the red solid lines that 
are overstriked with a white solid line show the 
PSF-affected input model. In panel b, the corre¬ 
sponding g — i colour measurements are shown 
with crossed circles, and the thick overstriked 
line shows the colour of the PSF-alfected model. 


region with a steeper decline (Freeman 1970). A third type was 
added to explain galaxies where surface-brightness profiles be¬ 
come shallower beyond a break at a large radius; these type-III 
profiles show an ‘anti-truncation’ of excess light that appears at 
3.2-6.0 inner scale lengths (Erwin et al. 2005). The inflection in 
type-III-s profiles - where the outer profile is believed to be part 
of the spheroid - is smooth and curved, and whilst isophotes are 
elliptical in the inner region for inclined galaxies, they become 
rounder at large radii. Meanwhile, the outer isophotes are not 
signficantly rounder in the other subclass, type Ill-d, where the 
outer profile is thought to be part of the disc. My selection of 
objects to smdy is made in Sect. 4.1, which also outlines the 
remaining parts of this section. 


4.1. Details of my selection, the observations, and the models 
of the five examined galaxies 

I selected five galaxies from EPB08, which all have outer parts 
that are classified as certain, uncertain, or possible type-III or 
type-III-s structures. Two galaxies are discussed in Sect. 4.2: 
NGC 7280 shows an odd structure and NGC 4102 was observed 
by two studies who report discrepant profiles and a possible type 
III-s structure. Three additional galaxies where the influence of 
scattered light is perhaps less clear are discussed in Appendix D.2: 
NGC 3507 and NGC 4477 were observed with SDSS and are 


classified with different types, and NGC 2880 is classified as type 
III-s. 

I modelled these galaxies using the same single-disc model 
parameters as EPB08.1 also added a bulge to each galaxy model, 
which I fitted by trial and error to have convolved surface- 
brightness structures match the (/?-band) measurements in the 
centre. Details of the observations are given in Table A.l and 
model parameters including galaxy types and inclination angles 
are collected In Table D.l. 

EPB08 discuss how the use of free or fixed fits of ellipses 
affect the profiles. They clearly show that free ellipse fits result in 
wrong profiles in the centre regions of the galaxies, as a free fit 
may track bar-distorted isophotes. I tilted all disc models using 
the provided inclinations to achieve elliptical isophotes, and there¬ 
after obtained the surface-brightness profiles with fixed ellipses. I 
never tilted the surface-brightness structure of the centre bulge. 

4.2. Two galaxies where scattered light dominates faint 
structures: NGC 7280 and NGC 4102 

Figure 4a shows surface-brightness profiles of my models of NGC 
7280 together with the /?-band measurements of EPB08. An anti¬ 
truncation break appears at rabr - 100". The convolved profiles 
overlap the measurements very closely out to the scattered-light 
halo radius 50 rno 60". Outside of both rno and rabr, there 
is more scattered light than the measurements indicate. The sug- 
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Fig. 4. Radial R-band and i-band surface- 
brightness profiles versus radius for two face- 
on disc galaxies: a) and b) NGC 7280, and c) 
and d) NGC 4102. Panels, lines, and coloured 
bullets are as described in Fig. 2. For NGC 
7280, the measured centre magnitudes at 1" 
seeing are 1.5 mag arcsec^^ fainter than the 
model magnitude at the centre; there is no ad¬ 
ditional offset here between the 1-band and the 
R-band profiles. Measurements of EPB08 are 
shown with crosses in panels a and c (for the R 
band, their fig. 14). The r-band measurements 
of Courteau (1996) for NGC 4102 are also 
shown with black dots in panel c. Red sym¬ 
bols indicate values that are fainter than the 
sky-uncertainty limit. 


gested type III structure seems to be explained by the procedure 
where EPB08 mask extended ‘haloes’ of nearby bright stars, 
whereby they also remove halo light from NGC 7280 (cf. their 
sect. 5). The colour profiles in Fig. 4b show red excess immedi¬ 
ately outside of the break radius at = 4275. For larger radii, 
the values are A(R - i) ^ 0.23 mag arcsec^^ at r = 60" and about 
0.84-0.96 mag arcsec^^ at the reported anti-truncation break ra¬ 
dius Tabr = 100". The amount of red excess decreases outside the 
maximum at r = 120". The increasingly red model profile for 
r ^ 135" originates in the outer parts of the bulge. 

Corresponding surface-brightness profiles of the models of 
NGC 4102 are shown in Fig. 4c. They show a good match be¬ 
tween the R-band measurements and the prohles of the convolved 
models that use either PSFyom or PSFkvi- The kink in the con¬ 
volved profiles of NGC 7280 and NGC 4102, where 30 r ^ 40", 
is due to the limiting radius where isophotes at larger radii were 
fitted with ellipses instead of circles. Using PSFkvi, the anti¬ 
truncation radius r^br ~ 100" agrees well with the scattered-light 
halo radius rno. 

The measurements of Courteau (1996) are brighter than those 
of EPB08 by up to 1 mag arcsec^^for r $ 150". A difference 
such as this could be explained by different PSFs or inaccurate 
sky subtraction. The figure suggests that the PSF of Courteau 
(1996) could be closer to PSFy^n,; his measurements are clearly 
underestimated for r ^ 150", where yu ^ 26.7 mag arcsec^^, but 
these data points are fainter than his reliability limit value /j. = 
26magarcsec“^. There is red excess immediately outside of the 
break radius at r^y - 4279 also for this galaxy, Fig. 4d. At larger 
radii, the values are A (R - /) 0.080-0.11 magarcsec^^ at r = 

60" and about 0.48 mag arcsec^^ at rabr- 

5. About analyses of galaxy structures that are 
fitted with Sersic-type profiles 

Some of the first astronomical studies that address the role of scat¬ 
tered light examine observations of nearby elliptical and lenticular 


galaxies that appear large, where these effects are found to be 
small and even insignificant. Galaxies included in such studies 
include NGC 221, NGC 3115, NGC 3379, NGC 4494, and NGC 
4649 (de Vaucouleurs 1948, 1953; de Vaucouleurs & Capaccioli 
1979; Capaccioli & de Vaucouleurs 1983; Capaccioli et al. 1987) 
and M 31 (de Vaucouleurs 1958). Excess light in the outer regions 
of the cD galaxy A 2029 closely follows de Vaucouleur’s law, 
after that the diffuse scattered-light component of field stars is 
carefully removed (Uson et al. 1991). A less dramatic influence 
of scattered light is found in a later study of envelopes of four 
additional cD and D galaxies (Mackie 1992, NGC 7647, NGC 
7720, NGC 7728, and A 407 Gl). 

Michard (1999) and Michard & Poulain (2000) correct colour 
prohles in centre regions of elliptical and lenticular galaxies from 
effects of variable PSFs (“differential seeing”). Idiart et al. (2002) 
study colour prohles of 36 elliptical and lenticular galaxies and 
correct colours for wavelength-dependent scattered-light effects; 
to remove artihcial colours, they convolve data of one hlter with 
the PSF of the other hlter. ^ M02 studies the lenticular galaxy NGC 
3115 as well as the four elliptical galaxies NGC 4406, NGC 4473, 
NGC 4551, and NGC 4550. The two latter galaxies are smaller 
than the two earlier elliptical galaxies, and are found to be more 
strongly affected by scattered light. The general opinion, however, 
that PSF effects are unimportant in elliptical-type galaxies, seems 
to originate in the early studies. 

The studies I mention above provide object-specihc results 
regarding the influence of scattered light. There is, so far, no 
general Sersic-proflle parameter study that shows how strong 
scattered-light effects are in such Sersic-type galaxies of variable 
properties. A systematic search of scattered-light haloes in a large 
parameter range that includes objects that appear both small and 
large is provided in Sect. 5.1. Host galaxies of BCGs are typically 
modelled with n = 1, and these objects are therefore relatively 


^ Studies that choose to use this method should check its validity as 
there is no comment about its accuracy for different sets of parameters. 
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strongly affected by scattered light. Multiple examples of BCGs 
are presented in Sect. 5.2. Thereafter, an example of a galaxy with 
intermediate effects of scattered light is discussed in Sect. 5.3, 
and two examples of large galaxies where effects are small are 
shown in Sect. 5.4. 


5.1. A systematic search for the appearance of scattered-light 
haloes in Sersic-profile galaxy models 

1 calculated and analysed a large number of elliptical-type galaxy 
models, to characterize PSF effects for a mostly complete range 
of the Sersic-prohle parameters n and r^. Eight values on n were 
used with the models: 0.5, 1, 1.5, 2, 3, 4, 6, and 8. For each 
value of n, a V-band and an /-band model were calculated using 
sixteen different half-light radii ry. 1, 2, 3, 4, 5, 7, 10, 15, 20, 25, 
30, 40, 50, 70, 100, and 200". Not all parameter combinations 
are expected in structures of real galaxies. Larger values on n 
and Te require grids with very large pixels, and therefore low 
resolution in the centre. I assumed that - p^ y — 2. Each 
model was convolved with PSEv/Qm^ PSEkvi, and PSF,om- The 
scattered-light halo radius ri m is shown versus and n for all 
three PSFs in Fig. 5. As the figure shows, a scattered-light halo 
appears at smaller radii and lower magnitudes for lower values of 
n and r^. 

With PSFypm (Fig. 5a) and when n ^ 2.0, haloes are visible at 
a magnitude that is <10 mag arcsec^^ fainter than the centre sur¬ 
face brightness, and a radius rno hr^. The corresponding radius 
is riio 4re (6.7 h) for galaxies where n = 1.0, which includes 
edge-on and face-on disc galaxies. Surface brightnesses at rno 
are for higher values on n and so faint that the scattered-light 
halo is unlikely to be observed in galaxies that are aheady faint 
at the centre; coarse limiting parameter values are, for example, 
n — 3.0 and ^ 3" ox n - 4.0 and re ^ 275. 

The scattered-light halo radius r\ lo moves outwards to larger 
radii with the lower-limit scattering PSFkvi (Fig. 5b); for example, 
riio ^ 4.5re for n = 1.0 and rno $ 9re for n — 2.0. Examples of 
coarse limiting parameters for models where surface brightnesses 
are too faint to be observed at rno are ^ 175 for n - 3.0 and 
Te ^ 1" for n = 4.0. In contrast, because of the red-halo effect 
(Sirianni et al. 1998), rno moves inwards to shorter radii with 
the more strongly scattering PSF, om (Fig. 5c). This is particularly 
noticable for values 10-20". 

5.2. On observations and general properties of BCGs 

Sargent & Searle (1970) conclude that compact galaxies are iso¬ 
lated extragalactic FIii regions. Later, Searle & Sargent (1972) 
and Searle et al. (1973) hnd two possibilities regarding the origin 
of BCGs; at first, they argue that they are young galaxies that 
form their hrst generation of stars, but they settle with an explana¬ 
tion where star formation occurs in short bursts that are separated 
by long quiescent periods. Various morphological classes are 
introduced to explain both the star-formation region and the en¬ 
closing host galaxy (Loose & Thuan 1986b; Kunth et al. 1988; 
Salzer et al. 1989; Telles et al. 1997). Radial isophotes of faint 
host galaxies are found to decline as an exponential (e.g. Bergvall 
1985; Telles & Terlevich 1997), a de Vaucoleur’s law (i.e. a Sersic 
profile with n - 4, e.g. Loose & Thuan 1986a; Doublier et al. 
1997), multi-component models (Papaderos et al. 1996, hereafter 
P96) and with Sersic prohles that use other values on n (Bergvall 
& Ostlin 2002, hereafter BO02; Caon et al. 2005, hereafter C05; 
Amorfn et al. 2007, 2009). General physical and observable prop¬ 


erties of BCGs are overviewed by, for example, Kunth & Ostlin 

( 2000 ). 

Comparisons of structural properties, average colours, and 
colour gradients between different classes of dwarf galaxies pro¬ 
vide a mean to test evolutionary scenarios between the classes 
(Thuan 1985; Bothun et al. 1986; Davies & Phillipps 1988). The 
assessment of properties of the host is important to establish the 
evolutionary state of the galaxy and its star-formation history. 
However, the host is fainter than the starburst regions, which 
makes it difficult to measure. 

In hosts with an exponential decline, scattered-light haloes ap¬ 
pear at relatively bright intensities and small radii ri lo - 6.1-1.5h 
(4-A.5rO, see Fig. 5. Integrated light from the host is modulated 
by scattered light from the bright, compact, and possibly asymmet¬ 
rically situated central starburst regions, which may complicate 
the interpretation of the structure. Scattered-light haloes should 
be visible in BCGs, but they were never reported. 

Here, I discuss surface-brightness prohles of hve BCGs, with 
a focus on the faint outer parts of the host. I selected galaxies 
where there are measurements for r ^ 25", to avoid problems 
with the relatively high seeing of the PSFs I use and resulting 
poorly resolved central regions. Most BCGs are very small and 
are only observed to shallow depths. I study surface-brightness 
prohles - including the role of the ellipticity - of the extreme BCG 
ESO 400-G043 in Sect. 5.2.1, which surface-brightness prohles 
are dominated by the bright starburst regions. C05 hnd that a 
Sersic prohle hts the measurements of Mrk 5 and I Zw 123 better 
than an exponential prohle. I discuss how scattered light could 
be a part of the explanation of the result for Mrk 5 in Sect. 5.2.2, 
which in comparison to ESO 400-G043 clearly also shows parts 
of the exponential prohle of the host. Three additional BCGs 
are discussed in Appendix E: the extreme BCG ESO 350-IG038 
(Haro 11) that hosts asymmetrically placed starburst regions; 
observations of the BCG Mrk 297, which seem to be affected by 
varying levels of sky subtraction in addition to scattered light; and 
UM 465, which is included in a recent extensive study of BCGs 
(Micheva et al. 2013b, hereafter Ml 3b). Parameters of observed 
and modelled surface-brightness prohles of the discussed BCGs 
are collected in Table E.l, and details of the observations are 
given in Table A. 1. 

5.2.1. The extreme BCG ESO 400-G043 

Out of four BCGs in BO02,1 htted the surface-brightness prohles 
of ESO 400-G043 with the smallest value on the scale length, 
resulting in effects of scattered light that are slightly larger than in 
ESO 350-IG038 and much larger than in ESO 338-IG04 or ESO 
480-IG12. BO02 present B-band prohles of ESO 400-G043 that 
they ht with both an exponential disc and a Sersic prohle. More 
recent V-band prohles are presented by Micheva et al. (2013 a, 
hereafter M13a), who measure an ellipticity of e = 0.22. The 
authors ht their measured surface-brightness prohles with two 
exponential discs in two separate magnitude ranges. 

I set the model parameters to achieve a reasonable agreement 
between convolved structures and both the B-band and the V-band 
measurements. My model consisted of a bulge, where n = 1.0 and 
Te = 175, that is superposed on an exponential-disc host where 
n - 1.0 and - 473, = 5.75 mag arcsec^^, and 

Be,B “ Pej - 0.5 mag arcsec^^. I also set the intensity constant for 
r < 5". Major-axis surface-brightness prohles, the measurements, 
and all three hts are shown in Fig. 6a. 

Outside a fairly hat core, intensities decrease steeply until 
they are about p = 5 mag arcsec^^ fainter than at the centre. A 
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Fig. 5. Scattered-light halo radius rno in Sersic-type galaxies that are modelled with the Sersic profile, rnof^ * is drawn versus for n = 0.5, 1.0, 1.5, 
2.0, 3.0, 4.0, 6.0, and 8.0. Each model was convolved using three different PSFs: a) PSEyom. b) PSEkti, and c) PSE,_om- The magnitude at rug relative 
to the maximum value 0 mag arcsec^^ at the centre is indicated with three different symbols: *,3 < yu < 7 mag arcsec^^; *,7 < /j < 11 mag arcsec^^; 
and o, 11 < ^ < 15 mag arcsec^^. The symbol sizes are scaled linearly, where larger (smaller) symbols indicate a lower (higher) magnitude. 


scattered-light halo dominates at fainter intensities, but begins 
already near the centre; see the locations of the coloured bullets 
in the figure. The Z?-band measurements match the profile of 
the model that was convolved using PSFyom well. The V-band 
measurements better match the profile that used PSFy 3 m. The 
dependence on scattered light with time - as is represented by 
the different PSFs - is clearly seen in the separation between 
the profiles of the convolved models. For example, compare the 
blue and the purple lines that differ by about 3 mag arcsec^^ at 
r = 32". My parameter values are lower than those of BO02 
and M13a; they measure their profiles in the radial range where 
scattered-light effects appear to dominate, r S 14". 

The B - / colour profiles in Fig. 6b clearly show more red 
excess with the two earlier PSFyom and PSF, om than with the 
later PSFv 3 ni and PSF, 3 ni. Examples of values on the red excess 
at r = 30" are A(B -1) ^ 1.6magarcsec“^ with PSFyom and 
PSFy^Oni, and about 1.0 mag arcsec ^ with PSFv; 3 m and PSF,^^. 

In view of the results of this study, the outer parts of the 
surface-brightness and colour profiles of the more distant ESO 
400-G043 and ESO 350-IG038 appear extreme because of the 
scattered light, which originates in the very bright starburst re¬ 
gions in the centre. The starburst regions in the two BCGs ESO 
338-IG04 and ESO 480-IG12 are fainter and the profiles thereby 
less extreme (not shown here). My study suggests that the host is 
visible in the inner regions of the extended parts of these two lat¬ 
ter BCGs as well as in Mrk 5 (see helow), whilst it is completely 
concealed in the first two objects. Effects of scattered light are 
significant in ah four galaxies due to the small value on re. I did 
not separately model the bands B, V, or K, but it seems reason¬ 
able that a similar conclusion applies to the remaining strange 
red-halo objects that are reported by Bergvall et al. (2005). 

Besides inducing haloes and red excess, scattered light affects 
the flattening that is measured in elongated objects. I show one 
example here to illustrate how measured ellipticities of objects 
such as ESO 400-G043 depend on the PSE. In this context, Ml3a 



r (arcsec) 


Fig. 7. Ellipticity versus radius e(r) of convolved models of ESO 400- 
G043; although, using the model ellipticity e = 0.45 (dash-tripple-dotted 
line). The model was convolved using five PSFs: PSFkti (purple line), 
PSFyon, and PSFysm (blue solid and dash-dotted lines), PSF,om and 
PSF,,3,„ (red solid and dash-dotted lines). 

measure the ellipticity e - 0.22 for this BCG. I instead used 
e = 0.45, to enhance the effects of scattered light. I convolved the 
flattened model with five PSFs and measured the ellipticity in the 
resulting model, see Fig. 7. 

None of the measured ellipticities reach the constant model 
value. All measurements show more circular values, i.e. the el¬ 
lipticity is smaller. PSFs with more scattered light show more 
circular structures; compare, for example, the curve of PSF/ ^m 
with that of PSFkvi. In the innermost region, the ellipticity in¬ 
creases with radius for all PSFs, because of the bulge that was 
not flattened. For example, using PSFyom, the ellipticity increases 
slightly from e(r ^ 3") = 0.25 to e(i5 ^ r ^ 20") = 0.37. At 
larger radii, the ellipticity decreases towards zero for all PSFs. 
Scattered-light effects in flattened galaxies are more pronounced 
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Fig. 6. Surface-brightness profiles versus radius 
for the BCG ESO 400-G043. Panels, lines, and 
coloured bullets are as described in Fig. 2. The 
/-band profiles are offset by -1.5 magarcsec^^ 
from the 6-band profiles for increased visibil¬ 
ity. Grey and black bullets and error bars mark 
measured 6-band (BO02) and V-band (M13a) 
values, respectively. The grey dotted (dashed) 
line shows the best-fit Sersic-type model (expo¬ 
nential fit to the disc in the range 16 < r < 32") 
of BO02. The black dotted line shows the best 
exponential fit to the disc of M13a (for 24 < 
Pv ^ 26magarcsec“^). 


on the minor axis than on the major axis, but measurements along 
ellipses deviate minutely from major-axis measurements (not 
shown). It appears that the only way to measure true radial ellip- 
ticity structures is to deconvolve the measured surface-brightness 
structure. 

5.2.2. The BCG Mrk 5 

Observations of Mrk 5 were made in the 6-band (Cairos et al. 
2001, hereafter COl), who fit the measurements with an expo¬ 
nential profile. The data are reassessed by C05, who fit the host- 
galaxy measurements with a Sersic profile; these authors extend 
the set of published data points to r = 58". Amorin et al. (2009) 
also use the data to make a two-dimensional Sersic-profile fit. 

I matched the measurements with a model using a bulge 
where n = 1, re = 2"0, and pe = 21.3 mag arcsec^^, and a host 
where n = 1, re = 8"0, and pe = 22.0 mag arcsec^^. I also set 
Pe.B -yUe,/ = 0.5 mag arcsec^^. It is only possibly to achieve a fair 
agreement between the bulge and the measurements for r 3", 
but there is only little intensity in this region, which effects on the 
larger object are small. Surface-brightness profiles of the models 
and the three sets of profile fits are shown in Fig. 8a. 

According to M13a, the PSFmb of the NOT is close to PSFkti 
and even decreases for r ^ 30", see Fig. A.l; they conclude it is 
too faint to influence the observations. No PSF shows a decline 
that is steeper than a r“^ power-law decline (also see figs. 1 and 2 


in Paper I) and it seems unlikely that the PSF of the NOT is any 
different. With my model parameters, the outer measurements 
of Mrk 5 appear to require a PSF that is as bright as PSFv^m to 
be explained as scattered light. There are three circumstances 
that make this possible. The NOT PSF could, as other PSFs, 
vary with time and acmally be as bright as PSFv 3 m, there is not 
enough information about which data were used to derive PSFmis- 
The PSF needs to be 110% as extended as the measurements to 
correctly estimate its influence (Appendix B), considering PSFmis 
this would affect all values for r ^ 25", including nearly all of 
the host (or halo). It is unknown how the used sky-subtraction 
procedure has affected the measurements, in particular where 
Pb > 26 mag arcsec^^ for r ^ 24" (compare with the discussion 
for Mrk 297 in Appendix E.2). 

The scattered-light halo radius rno - 25" for PSFyom and 
rno - 36" for PSFkti, whilst it is closer to rno - 22" for 
PSFi/ 3 in. The 6-band measurements of COl and the reassessed 
values of C05 overlap each other and the convolved model using 
PSFy 3 n,; the values deviate slightly for r ^ 3", but that deviation 
is unimportant to the discussion here. The exponential fit of COl 
agrees with the value I chose for the single-disc host component. 
C05 fit values at all radii with a single Sersic profile and thereby 
need to use the higher value n - 2.83; they also use an equally 
high value with the smaller appearing BCG I Zw 123 (not shown 
here). The fit of Amorfn et al. (2009) agrees poorly with the 
measurements in large intervals (as with Mrk 297). 
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Fig. 8. Surface-brightness profiles versus radius 
for the BCG Mrk 5. Panels, lines, and coloured 
bullets are as described in Fig. 2. The /-band pro¬ 
files are offset by -l.Smagarcsec^^ from the 
S-band profiles for increased visibility. Black 
(grey) diamonds and error bars mark measured 
B-band values of COl (C05). The black (grey) 
dotted line shows the best-fit Sersic-type model 
of COl (COS); the black line is difficult to see, as 
it falls atop the /J-band model line. The dash- 
tripple-dotted line shows the fitted profile of 
Amon'n et al. (2009). 


The B - I colour profiles in Fig. 8b show dramatic amounts 
of red excess for r ^ 18". For example, at r = 30" the red excess 
is A (B - /) ^ 0.31-0.37 mag arcsec^^ and at r = 40" the values 
are A(B - /) ^ 0.8-1.0 mag arcsec“^. 


5.3. About galaxies where scattered-light effects are 
intermediate 

I classify intermediate effects of scattered light as those that ap¬ 
pear in surface-brightness stmctures where the brighter parts are 
well fitted with 1 < « < 4.1 study one such galaxy here, NGC 
4551. Kormendy et al. (2009, hereafter K09) present observations 
of NGC 4551 of eight different sources in the B, V, and R bands. 
The authors fit a V-band model to the data with the parameters 
n = 1.98, Te = 1575, yUe = 20.75 mag arcsec^^, and D - 20Mpc. 
The profile extends out to r = 110" where p ^ 27.3 mag arcsec^^. 
M02 presents observations that extend out to r = 70". 1 modelled 
NGC 4551 using n - 2.0 and re = 15", neglecting the somewhat 
brighter region where r < 278. To broaden the applicability of 
this one study to faint structures of similar but brighter galaxies, 
I extended the models out to r = 300". Surface-brightness and 
colour profiles as well as the V-band measurements are shown in 
Fig. 9a. Other Virgo-cluster galaxies, at about the same distance, 
which surface-brightness profiles are fitted with similar param¬ 
eters, include IC 0798, IC 3509, NGC 4387, NGC 4434, NGC 


4458, NGC 4478, NGC 4515, VCC 1407, VCC 1545, VCC 1828, 
and VCC 1910 (K09). 

The scattered-light halo radii are rno - 80" (PSTyom); ~ 
105" (PSFkvi), and ~ 45" (PSF,,om), also see Fig. 5. The first two 
radii are larger than (similar to) the radial range that is covered 
by M02 (K09). Differences between profiles that were calculated 
using PSFk7i, PSFyom, and PSFy 3 m (or PSF,,om and PSF, 3 m) are 
small, and are more pronounced where r > 80". Considering 
that the PSFs are less accurate at large radii, say for r ^ 100", 
it is difficult to draw a conclusion on how sensitive the outer 
structure is to temporal variations in the PSFs. Convolved surface- 
brightness profiles of the separate PSFs differ at r = 80" by 
at most 0.15 mag arcsec^^ in both bands, and at r = 150" by 
0.3 mag arcsec^^ in the V band and by 0.5 mag arcsec^^ in 
the I band. The convolved V - I colour gradient differs from the 
model gradient throughout the radial range, also where r < rno, 
cf. Fig. 9b. There is red excess for r > 25", which reaches 
V-I ^ 0.29magarcsec^^ at r = 130" using PSFyom and PSF/om, 
and V - I ^ 0.30 mag arcsec'^ using PSFy 3 in and PSF, 3 m. The 
colour gradients are positive in the regions 20 ^ r ^ 150" and 
r ^ 250". The noticeable bump between the two colour profiles, 
where 150 ^ r ^ 250", appears as a result of how I extrapolated 
the PSFs at larger radii. 

Surface-brightness profiles of galaxies fitted with models 
using similar Sersic parameters, such as those that I mention 
above, will show similar amounts of red excess. Whilst Scattered- 
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Fig. 9. Surface-brightness profiles verus radius for NGC 4551. Panels, 
lines, coloured bullets, and offsets are as described in Figs. 2 and 4. The 
/-band profile in panel a is offset by -0.5 mag arcsec^^ from the V-band 
profile. Black bullets show a fitted Sersic profile (the parameters are 
taken from K09). 

light effects always appear at some radius, the measured surface 
brightness is often too faint to reveal the scattered-light halo. 
However, the red excess seen in colours is detectable also at 
lower magnitudes and radii. Galaxy profiles fitted with smaller 
have steeper colour gradients and brighter scattered-light haloes; 
examples of such galaxies in the Virgo cluster include IC 0809, 
IC 3461, IC 3470, IC 3490, IC 3653, NGC 4464, NGC 4486A, 
NGC 4486B, VCC 1871, NGC 4467, VCC 1199, VCC 1440, and 
VCC 1627 (K09). 

5.4. About galaxies where scattered-light effects are smaller 

Scattered-light effects in surface-brightness structures of larger 
elliptical and cD galaxies (« > 4) are typically small to very small. 
Instead of discussing the appearance of scattered-light haloes at 
large radii and faint intensities in such galaxies, I focus on colour 
gradients at smaller radii where r < 250". Surface-brightness 
profiles are shown in Figs. 10a and 10c for two larger galaxies, 
NGC 4473 and NGC 4374. 

K09 fit observed surface-brightness structures of NGC 4473 
with n = 4.0, r^ = 5178, using the distance D - 15.28 Mpc, and 
of NGC 4374 with n = 8.0, = 142", and D = 18.45 Mpc. I 

used the same model parameters, and I set r^j — 0.91 r^ y for 
both galaxies. With these galaxy models, a scattered-light halo 
appears at ^ 600" and rf°'^4374 ^ 3200 " (cf. Fig. 5 

and Sect. 5.1). Both convolved models show V -1 colour profiles 
that deviate from the respective model profile where r ^ 250", cf. 
Figs. 10b and lOd. 

The radial colour profile of NGC 4473 increases with radius 
for r ^ 80" - in agreement with the result of M02 for the same 
galaxy - and decreases for larger radii. Intensity differences due 


to the used temporally variable PSFs are minute; colour profiles of 
the convolved models differ by < 0.01 magarcsec^^ for r < 80" 
and by < 0.28 mag arcsec^^ for r < 250". Colours are much less 
affected than with lower Sersic-profile parameter values, but there 
is still some blue excess in the innermost region and red excess 
in the remaining region. The maximum colour difference across 
the shown radial range is (V - /)max - 0-08 mag arcsec^^. The 
agreement is striking between these colour profiles and those of 
the much larger galaxy NGC 4374, both regarding amplitudes 
and shape. It seems that in this case the colour difference appears 
because of the convolution with a PSF, but the result is nearly 
insensitive to the exact properties of the model. 

Compared to models of smaller galaxies with low values 
on n, it is comparatively easy to deconvolve measured data of 
large galaxies such as these. The accuracy requirements on the 
measurements and the PSF to remove the colour artefacts are 
small. 


6. Discussion 

The presence of diffuse excess light around galaxies is ubiquitous, 
and it is tightly linked with the existence of colour excess at 
intermediate to large radii. The outcome and implications of this 
paper are discussed in the following sections. 

6.1. Edge-on galaxies 

All three scrutinies of haloes around edge-on galaxies studied 
here, including IC 5249, NGC 4565, and NGC 4244, corroborate 
the result for NGC 5907 in Paper I - the faint parts of the ob¬ 
served surface-brightness profiles could be explained by scattered 
light. Specifically, Fry et al. (1999, hereafter F99) fit their surface- 
brightness profiles of NGC 4244 with a single exponential disc 
and argue that the galaxy does not have a halo. My scrutiny shows 
that their data trace the expected scattered-light profiles at large 
radii, where they deviate from an exponential profile. The agree¬ 
ment between photometry and resolved-stars measurements of 
this object is poor, which illustrates the need for further research 
to understand the differences. 

A number of scenarios were proposed to explain the forma¬ 
tion of currently discovered thick discs; see the overview in, for 
example, Yoachim & Dalcanton (2008). It is argued that thick 
discs are a necessary consequence of disc formation in galaxy 
formation models (Comeron et al. 2011), and they are described 
as a ‘local missing baryons’ reservoir. Another study discusses 
the thin-disc to thick-disc dichotomy in the Milky Way, and ar¬ 
gues for the term ‘thicker disc component’ instead of ‘thick disc’ 
(Bovy et al. 2012). Both the thick disc and the halo indicate a 
presence of excess light, and they can be treated similarly (Morri¬ 
son 1999). My results regarding diffuse thick discs is no different 
to the haloes. The two thick discs of FGC 310 and FGC 1285 
bracket all the 34 objects of YD06 (the effects of scattered light 
are stronger in FGC 1063 than in FGC 310 since its apparent 
size is smaller). All evidence points at scattered light being the 
ubiquitous component that conceals any real thick disc in addition 
to modifying the scale height of the thin disc. The geometrical 
shape of surface-brightness profiles of edge-on galaxies implies 
strong effects of scattered light in faint regions, regardless of the 
size and scale height of the studied object. 

The analysis of the role of scattered light in stacked data can 
be more difficult when exact dates of the galaxy observations are 
unknown and their geometrical properties differ. I analysed parts 
of the stacked LSBG data of BZCIO, who deconvolve their data. 
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Fig. 10. Scattered-light effects of two Virgo- 
cluster galaxy models: NGC 4473 (left panels) 
and NGC 4374 (right panels). Panels, lines, and 
offsets are as described in Figs. 2 and 4. The con¬ 
volved profiles are invisible as they fall atop the 
model profiles. Black bullets show a fitted Sersic 
profile for the respective galaxy (the parameters 
are taken from K09). 


but neglect temporally varying PSFs. My results show that it is 
possible to model colour profiles with similar amounts of red ex¬ 
cess in the outer parts as the presented deconvolved observations. 

Feldmeier et al. (2013) stack images of Lya emitting galaxies 
at z ~ 2.1 (data from Guaita et al. 2010) and at z ~ 3.1 (data from 
Gronwall et al. 2007; Matsuda et al. 2012). They extract PSFs 
from bright stars in their stacks and use them to deconvolve the 
respective stack. The objects are very distant and therefore appear 
as point sources, which is why compact PSFs suffice. The z ~ 2.1 
data were observed within eleven days (2007 December 3-13), 
which means that one average PSF should suffice to analyse the 
data accurately. The data at z ~ 3.1 were observed during a time 
interval of about 16 months (Gronwall et al. 2007) and 38 months 
(Matsuda et al. 2011, 2012), which means that a single PSF might 
be an inaccurate approximation. Whilst Feldmeier et al. (2013) 
explore an alternative origin of the halo for this case in other 
analysis-related issues, it cannot be excluded that the PSFs of 
the images used to find a (weakened) halo were slightly different 
from the mean PSF; the PSFs presented here in Fig. A.l show 
that differences are very small at the short radii of these point 
sources. 

6.2. Face-on disc galaxies 

Both PSFyom and PSFkvi appear to be close to the true PSF in 
the observations of NGC 4102 of EPB08. The surface-brightness 
profiles of the three galaxies observed with the SDSS and the 
WIYN telescope (NGC 4477, NGC 3507, and NGC 2880) also 
appear to be explainable with scattered light, except the outer¬ 
most faintest measurements. In each case, the scattered-light halo 
radius ri lo agrees well with the anti-truncation radius rabr, which 
is used as an indicator for type III-s profiles. It is important to 
remember that the models are simplified representations of the 
surface-brightness profiles of the real galaxies. Real variations 
in the intensity due to dust lanes or any azimuthal irregulari¬ 
ties are neglected. It cannot be excluded that PSFy^m is closer 


to the true PSF, and that my models are lower limits. Besides 
the too large subtracted sky component, the differences of up to 
about 1.5 mag arcsec^^ in the surface-brightness profiles of NGC 
7280 between the profiles convolved using PSFi/om, PSFv^m, and 
PSFk 7 i, indicate that the structures of this smaller appearing 
galaxy are highly sensitive to temporally variable PSFs. 

There are three further complications with the observed 
galaxy profiles: the PSFs are unknown, especially at larger radii, 
which makes it impossible to accurately dismiss their influence; 
the accuracy of the sky subtraction seems to differ in some cases; 
and the adopted eifipticity affects the resufting profifes and the 
anafysis. For exampfe, scattered light causes a decrease of the 
ellipticity at larger radii in a galaxy model with a radially con¬ 
stant ellipticity (see the study of the BGC ESO 400-G043 in 
Sect. 5.2.1). Inclination angles may be underestimated when data 
are not first deconvolved. 

The models of all five galaxies show significant amounts of 
red excess, which only appears outside of the inner break in those 
galaxies that have one. For galaxies without a break radius, the 
red excess begins already near the centre. Any study of radial 
colour gradients in the outer fainter parts of face-on disc galaxies 
is affected by scattered light. 

Not all observations of face-on disc galaxies are deep enough 
to detect a type III profile, and depending on the model parameters 
of the inner and the outer disc the anti truncation may only appear 
at very large radii and faint intensities. The models presented 
here present good evidence that excess light at large radii around 
face-on disc galaxies is scattered light; the type III profile and in 
particular the type III-s profile could be a sign for it. Only after 
the scattered light is removed, it is possible to see if the excess 
light instead comes from the extended parts of the bulge in the 
centre, as is currently assumed for type III-s. 

There are few identifications of haloes and anomalous colours 
in space-based observations. The uncertainties in the outermost 
parts of the HST PSFs are undetermined, this concerns both the 
Hubble Ultra Deep Field (HUDF; Beckwith et al. 2006) edge-on 
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galaxy study of Zibetti & Ferguson (2004) that is addressed by 
JOS, and the analysis of the haloes around UDF 3372 and UDF 
5417 by Trujillo & Bakos (2013). Observations with space tele¬ 
scopes, such as the HST, are less affected by temporally varying 
PSFs than ground-based observations, but it cannot be excluded 
that also here may there be some excess light at larger radii of the 
PSFs, so far, unaccounted for (cf. appendix A in Paper I). 

6.3. Galaxies fitted with Sersic-type profiles 

All galaxy profiles modelled with a Sersic profile are always - 
beginning at some radius - affected by scattered light, for all rea¬ 
sonable values on n and rg. The exact values on the radius and the 
surface brightness where the scattered-light halo appears depend 
on the object properties and the PSF. With higher parameter val¬ 
ues, the halo moves outwards to larger radii and fainter intensities, 
and resulting colour gradients become shallower. Whilst haloes 
of the largest appearing galaxies (n ^ 4) are so faint that they 
are unlikely to be observed, colour gradients are still somewhat 
affected in the centre region, as is seen here with NGC 4473 
and NGC 4374. At the other end of smaller appearing galaxies, 
a scattered-light halo appears at measurable magnitudes for all 
values of r^. when « = 1, and colours are dominated by differences 
between PSFs in large radial intervals. 

Numerous measurements of BCG hosts are affected by scat¬ 
tered light. My analyses of the extreme BCGs ESO 400-G043 and 
ESO 350-1G038 indicate that the host-galaxy part in the observa¬ 
tions is completely dominated by scattered light that originates in 
the bright starburst regions in the centre. The BCGs Mrk 297 and 
UM 465 were observed in different studies, they contain fainter 
starbursts and are less affected; here, the host galaxy appears to 
be visible near the centre, whilst scattered light still dominates in 
the fainter parts of the prohles. Current measurements of Mrk 297 
differ in the outermost parts, plausibly so because of differences 
in the sky subtraction or the PSEs at the time of the observations - 
results are ambiguous when both the sky subtraction is uncertain 
and the PSEs are unknown. The analysis of the BCG Mrk 5 that 
C05 present shows values of the Sersic index n ^ 2.6. The authors 
argue that this BCG together with 1 Zw 123 represents a new class 
of objects where indices n > 1. My analysis suggests that a single 
exponential profile suffices to explain the measurements when 
scattered light is considered. 

6.4. Effects of scattered light in other wavelength ranges 

This study has focused on effects of scattered light in the visual 
wavelength range. It has also earlier been found and studied at 
radio wavelengths, where it is referred to as ‘stray radiation’ (for 
example, Kalberla et al. 1980). 

Several studies present surface-brightness profiles for the 
same types of galaxies strudied here in infrared bandpasses (see 
the respective overview sections). As is also noted in Paper I, the 
PSFs of these bands, as measured with data of the 2MASS Large 
Galaxies Atlas (Jarrett et al. 2003, where observations were made 
at Mt. Hopkins in Arizona and the Cerro Tololo Inter-American 
Observatory), are even brighter at larger radii than the /-band PSE 
(Michard 2007). Extended PSEs are even less well determined in 
the infrared. Haloes, thick discs, and colour gradients are thereby 
just as questionable as those in the visual wavelength range. 

A recent study reports on the ubiquity of diffuse haloes around 
galaxies in the ultraviolet wavelength range (Hodges-Kluck & 
Bregman 2014). The study uses the Galaxy Evolution Explorer 
50cm telescope and the space-based Swift Observatory. The au- 
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Fig. 11. Radial surface-brightness profiles. Panel a) shows the model 
of Zheng et al. (2014) convolved with PSFzis.g (blue lines) and with 
PSFk 7 i (purple lines), respectively. The regular (one-dimensional) model 
convolved with the regular (one-dimensional) PSF is shown with a solid 
(dash-dotted) line. Panel b) shows the difference between the convolved 
regular and one-dimensional model. 

thors subtract some scattered light, but neglect integrated light 
altogether. Considering the results that are presented here in the 
visual wavelength range, it would seem obvious that scattered 
light plays an important role to explain the ubiquity also in the 
ultraviolet. 

6.5. Effects of scattered light are underestimated in analyses 
made in one dimension 

A recent study finds that the outer discs and haloes in its obser¬ 
vations of 698 disc galaxies are little affected by scattered light 
(Zheng et al. 2014). Their PSFs appear similar to PSFk 7 i (their 
figs. 28 and 29). The region they study extends out to about 10" 
(Zheng, priv. comm 2014). The authors present a one dimensional 
galaxy model convolved with their one dimensional PSFs that is 
very similar to the input model, and they conclude that scattered 
light has minor effects on their results. 

I modelled the disc galaxy in the g band both in one and two 
dimensions, using = 3'.'5 and a bulge where n = 1, rg = 075, 
and po-pe -4.0 mag arcsec^^; I matched their convolved model 
with my one-dimensional convolved model. The results depend 
on if the model is convolved with their PSF^is ^ or with, for 
example, PSFkti, in particular at the seeing-dependent centre, 
see Fig. 11. Differences between the one and two-dimensional 
treatments are larger with the better-seeing PSFzis,^ than with 
PSFk 7 i. At larger radii, r ^ 13", the comparison with the profiles 
of the lower-limit PSFk 7 i suggests that their gr-band PSF underes¬ 
timates the scattered light. Notably, the difference between a two- 
dimensional and one-dimensional treatment increases at larger 
radii for larger galaxies (compare the lines of the one-dimensional 
example models in Fig. A.2 with the other models); the reason 
is that more points contribute to the intensity at any one point in 
a two-dimensional approach compared with a one-dimensional 
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line. The amount of scattered light is in general underestimated 
at larger radii when the analysis is made in one dimension, and 
the analysis therefore needs to be made in two dimensions. 

6.6. Suggestions for future observations of surface-brightness 
profiles that, in particular, include haloes and thick discs 

My analysis has shown that we know little if anything about 
haloes and thick discs observed in diffuse emission. Time- and 
instrument-specific scattered-light effects need to be removed 
at first, before the existence of any real (red) halo is studied, 
measurements of host-galaxies are reassessed, or properties of 
stellar populations and amounts of non-baryonic dark matter and 
other phenomena are scrutinized (for example, BO02; Zackrisson 
et al. 2006; Zackrisson & Flynn 2008; Zackrisson et al. 2012; 
Hernandez et al. 2013, M13a; M13b; Martm-Navarro et al. 2014). 
It also appears premature to discuss relations between thin-disc 
and thick-disc properties in edge-on disc galaxies (for example, 
Dalcanton & Bernstein 2002; Yoachim & Dalcanton 2005) be¬ 
fore the observations are corrected for scattered light; the same 
conclusion applies to kinematic studies of thick discs (Yoachim 
& Dalcanton 2005, 2008). 

The radius and intensity where haloes and thick discs appear 
depend on the object properties, the brightness of the sky, and 
the instrument and the telescope (as well as the atmosphere) 
that give rise to the PSF. The following points should be taken 
into account in future observations of haloes and thick discs to 
correctly address effects of scattered light: 

- The level of accuracy of the measurements must be deter¬ 
mined for each image or data set individually. Images that 
need to be stacked could be sorted in groups that span a time 
interval that is short enough to ensure that the PSF is invariant 
- say, two weeks.^ It is particularly important to determine the 
accuracy of the sky background if parts of the object structure 
are fainter than the sky. The accuracy can be tested by adding 
an exponential profile in different parts of the raw image, 
which is then reduced and analysed to find the level where 
the profile deviates from the exponential profile. My example 
models indicate that the accuracy of the background level 
of the sky is about 2.5 magarcsec^^ higher than the surface 
brightness at the radius where the reduced profile deviates by 
more than b = 0.10 mag arcsec^^ from an exponential profile. 

- In case of Sersic-type and face-on disc (edge-on) galaxy pro¬ 
files, it is necessary to use PSFs that are accurately determined 
out to radii that are 1.1 (1.5) times the maximum considered 
radius (vertical distance). It must be remembered that the PSF 
varies with both the wavelength and the time. Guidelines on 
how the PSF can be measured accurately out to large radii are 
provided in Paper I. 

- A model of the galaxy surface-brightness structure must ac¬ 
count for both the accuracy of the sky background and the 
contribution of scattered light, which define a maximum use¬ 
ful radius. Measurements must be brighter than the level 
where a convolved profile deviates from the modelled expo¬ 
nential profile by less than 6 - 0.05-0.15 magarcsec^^, and 
they must be inside of the maximum useful radius. Nothing 
can be said about measurements at larger radii. The analysis 
needs to be done more accurately if the measurements and 
the scattered-light model are similar and the intention is to 
use the data for scientific conclusions. 

I suggest a value this small to be on the safe side since the period and 
amplitude of actual temporal variations are unknown. 


- Points for further improvement of the analysis include ac¬ 
counting for possible flaws of the model, such as asymmetries 
and other data-reduction issues. It may also be necessary to 
use a two-dimensional PSF that accounts for spatial variations 
across the observed field. 

The ultimate solution is to deconvolve the data, but this ap¬ 
proach requires that both the measurements and the PSF are 
sufficiently accurate; it still remains to be determined exactly how 
accurate the data must be to do this. The required accuracy of the 
measurements, the sky level, and the PSF becomes prohibitively 
high already at intermediate radii for many galaxies. It is then 
difficult and even impossible to show that there is anything but 
scattered light at those radii. Instead, it appears more fruitful to 
analyse the data at smaller radii and then progress outwards as 
far as the errors allow. Except features such as much fainter tidal 
streams and tails, there should be little excess at larger radii. 

7. Conclusions 

I have presented multiple examples of models and measurements 
of surface-brightness structures of different types of galaxies. 
The results show that scattered light produces structures that 
are so similar to the faint measurements at large radii that they 
can explain their origin. Specifically, the affected parts of the 
structures include haloes and thick discs around edge-on galaxies, 
haloes and type III profiles around face-on disc galaxies, halo¬ 
like hosts around BCGs, haloes around elliptical galaxies, and 
increased scale lengths in the centre regions for smaller appearing 
galaxies of all kinds. 

The influence of scattered light has been modelled in terms 
of the minimum and larger amounts. As in Paper I, PSFk 7 i (King 
1971) has been used as a lower limit and the four brighter PSFyom, 
PSFy 3 ni, PSFj Om, and PSF, 3 n, of M02 as intermediate-to-higher 
amount PSFs - in particular. Paper I shows that PSFk 7 i can be 
used as a generally applicable lower limit with all instruments and 
(reflective) telescopes. With these models, it has been possible to 
study effects of temporally variable scattered light using radially 
extended PSFs in two bandpasses and one colour. All differences 
have been neglected between PSFs in the visual wavelength range, 
except in comparison with the I band. The surface-brightness 
structures have been modelled using standard galaxy-specific 
profiles. Each of the models has been convolved individually 
with all five PSFs, and the resulting surface-brightness profiles 
have been compared with measurements of different studies in 
different bandpasses. The combined influence of scattered light 
and the accuracy of the sky subtraction has also been analysed. 
Together, these two components provide strict limits on how far 
out it is possible to measure surface-brightness structures that can 
be accurately corrected. 

The colour profiles of all modelled galaxies show significant 
to drastic amounts of red excess across large radial intervals re¬ 
gardless of the galaxy type, and in particular in the faint halo 
regions. I have only studied one colour, such as R- i, but the re¬ 
sults are indicative of PSF-dependent colours in the outer regions. 
Any radial differences between PSFs of different wavelengths or 
bandpasses will induce artificial colour profiles that could be ei¬ 
ther red or blue. Notably, artificial colours are also seen at smaller 
radii than where scattered-light haloes appear. 

The selected sample of galaxy observations demonstrates 
that scattered light plays an integral part to explain faint surface- 
brightness and colour profiles around various types of galax¬ 
ies. Currently, diffuse data only show ubiquitous scattered-light 
haloes. To measure physical haloes it is necessary to first remove 
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the scattered-light component, which requires accurate measure¬ 
ments of extended PSFs (cf. Paper I) and accurate analysis as 
shown here. For as long as scattered light is not excluded as the 
cause of the excess light, it cannot be expected that the faintest 
parts of the modelled profiles are anything but scattered light, as 
it is a more likely solution than other more exotic explanations. 
Notably, the increasing number of deep surveys have so far rarely 
measured extended PSFs. There are hints that PSFs of these sur¬ 
veys may be brighter than the ones used in the analysis here (Due 
et al. 2015). Such bright PSFs would further amplify the effects 
and make the results even more alarming. Meanwhile, it seems 
that the only way to measure haloes are through observations of 
asymmetric tidal streams and tails as well as of resolved stars of 
nearby galaxies. I have limited this study to galaxies in the vi¬ 
sual wavelength range, but the result applies to other wavelength 
ranges, haloes, excess light, and colour profiles around all kinds 
of objects on the sky. 
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Table A.2. Example model parameters. 


Galaxy model 

n, Te / hi, ho, Apo 

'■tr.l 

rtr,2 

Fig. 

small Sersic 

1, 1.5 

4.0 

10.0 

A.2a 

medium Sersic 

1, 4.0 

10 

25.0 

A.2b 

large Sersic 

1,40 

10 

250 

A.2c 

small face-on 

13.5, 0.897,3.0 

4.0 

10.0 

A.2d 

medium face-on 

35.9, 2.39, 3.0 

10 

25.0 

A.2e 

large face-on 

359, 23.9, 3.0 

10 

250 

A.2f 


Notes. Column 1, example galaxy model; Col. 2, Sersic-profile index 
and half-light radius (arcsec), or the face-on galaxy inner and outer 
scale lengths (arcsec), and the intensity difference A//o = f^o.i - yt^o.o 
(mag arcsec^^); Col. 3, first truncation radius (arcsec); Col. 4, second 
truncation radius (arcsec); and Col. 5, results are shown in this figure. 
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Appendix A: Complementary measurements of the 
PSF and details of the observations 

The PSFs discussed in this paper are shown in Fig. A.l that 
is complementary to figs. 1 and 2 in Paper I. The plot includes 
PSFk 71, PSFyom, PSFi/ym, PSFi_ 0 m, PSF; 3 m, PSFmbh (Morrisou 
et al. 1994), and PSF“gjj (Paper I) that are all introduced in 
Paper!. The following three PSFs are studied here: PSFpgg (F99), 
PSFa 99 (Abe et al. 1999, hereafter ABC99), and PSFmb (M13a). 

Details of all observations related to the measurements scru¬ 
tinized in this paper are collected in Table A.l; the information 
allows a straightforward comparison of the measurements in the 
context of scattered light. The table also contains the details of 
the observations of NGC 5907 discussed in Paper I. 

Appendix B: Application of the method on example 
models of Sersic-type and face-on disc galaxies 

The example models of surface-brightness profiles of edge-on 
galaxies in Paper I are complemented here with models of Sersic- 
type and face-on disc galaxies. I delimit the former type galaxy 
examples to n = 1, which results in exponential slopes; galaxy 
profiles of this type are less affected when n » 1 (cf. Sect. 5). 
Galaxy stmctures well represented by Sersic-type profiles are less 
affected by scattered light than galaxies viewed edge on; simply 
viewed, the brightest intensities are concentrated to a single point 
instead of a line. Surface-brightness profiles of face-on disc galax¬ 
ies are also modelled with exponential slopes. However, effects 
of scattered light in the outer regions can be enhanced compared 
to my Sersic-type profiles when intensity declines in the central 
region are shallower inside a break radius; due to, for example, 
spiral arms or a bar. 


I calculated three example models of a small, an intermediate, 
and a large object for both the Sersic-type galaxy and the face-on 
disc galaxy. I varied the half-light radius and set the Sersic 
index n - I in the Sersic-type galaxy models. In the face-on 
disc galaxy models, I varied the outer scale length ho, set the 
inner scale length hi = l5ho, po,i = Afo.o + 3 mag arcsec^^, and the 
break radius rbr = 3ho', I did not model an additional bulge. Each 
model was convolved with PSFk 7 i and a radially extrapolated 
PSFyom - where both PSFs are more than two times as extended 
as the respective galaxy model. Additionally, each model was 
convolved with either PSF, after the PSF was truncated at the 
smaller radius i or the larger radius rtr_ 2 - Model parameters 
and truncation radii are given for each of the six models in 
Table A.2. The resulting surface-brightness profiles are shown 
in Fig. A.2. (Models were also convolved in one dimension with 
the respective PSF, the corresponding profiles are discussed in 
Secf. 6.6.) Values quoted in parenfheses below were modeled 
with PSFk 7 i, and the other values were modeled with PSFyom-1 
describe each of the three Sersic-type galaxy models separately, 
where I compare the resulting properties with the results of the 
respective face-on and edge-on disc galaxy model. 


Appendix B. 1: The two small galaxies 

The convolved profiles of the small galaxy show excess light at 
all distances r ^ 2'.'0, Fig. A.2a. Less light is seen within this 
radius; this decrease is largely caused by the centre part of the 
poor-seeing PSFs I used (FWHM ^ 2-A"). The scattered-light 
halo of the convolved model is too faint for r ^ 477 (r ^ 573) 
when the PSF is truncated at rtr = 4". The corresponding value 
for rtr = 10" is r ^ 973 (r ^ 973). The surface-brightnesses that 
were calculated using the three versions of the truncated PSFyom 
differ by up to abouf 0.7 magarcsec^^ af fhe cenfre, due to the 
variations of the individual normalization of the PSFs. The cor¬ 
rect amount of scattered light at r = 15" is only achieved when 
the PSF is not truncated within, say, r ss 16" (15" x 10"/973). 
Compared to the centre, both PSFk 7 i and PSFyom are about 
lOmagarcsec^^ fainter at r = 15" (see Fig. A.l). The model 
surface brightness is about 4.7 (3.3) magarcsec“^ fainter than 
the convolved stmcture at r = 12", where the surface brightness 
is Ap - 14 mag arcsec^^ fainter than at the centre; the signal-to- 
noise (S/N) value that would be required to extract the model 
structure from the observed (convolved) structure is about 76 
(21). At r = 775 (Ap ^ 8.7 mag arcsec^^), the corresponding 
values are 2.0 mag arcsec^^ and S/Nss6.2 (0.95 mag arcsec^^ and 
S/N=s2.4). 

The face-on galaxy shows excess light for r ^ 377, which is 
just outside of the break radius rbr = 277, Fig. A.2d. The centre 
plateau results in brighter profiles al larger radii. The model 
surface brighlness is abouf 5.8 (4.3) magarcsec^^ fainter lhan 
fhe convolved slruclure al r = 14" (Ap ^ 14magarcsec“^); fhe 
S/N value lhal would be required to exlracl fhe model struclure 
from the observed structure is about 200 (53). At r - 977 (Ap ^ 
8.7 mag arcsec^^), the corresponding values are 2.4 mag arcsec^^ 
and S/Nss9.2 (1.2magarcsec“^ and S/Nss3.0). These values are 
0.25-1.1 mag arcsec^^ brighter than for the Sersic-type galaxy. 

For the small edge-on galaxy model in Paper I, the model 
surface brightness is about 6.7 (5.1) magarcsec^^ fainter than 
the convolved structure atz = 1077 (Ap ^ 14magarcsec“^). At 
z = 771 (Ap ^ 8.7 mag arcsec^^), the corresponding value is 3.3 
(1.9) magarcsec^^. Compared to the Sersic-type galaxy, these 
values are 0.95-2.0 mag arcsec^^ brighter. 
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Fig. A.l. PSF surface-brightness profiles ver¬ 
sus the radius r for a Omag point source. In¬ 
dividual PSFs are drawn with coloured lines 
and symbols as indicated in the figure. PSFkti, 
PSFi/om. PSF, On,, PSF, 3 n,, PSFmbh. and PSF^g^ 
are shown in this figure as references. The light 
grey region indicates the PSF core, and the white 
region the aureole, as defined in Paper I. The ex¬ 
trapolated part of PSFa 99 is shown with a dotted 
line. 


Appendix B.2: The two intermediate-size gaiaxies 

There appears to be less excess light for the intermediate-size 
galaxy, at least near the centre, Fig. A.2b. The convolved profiles 
show excess light at all radii r ^ 5". When the PSF is trun¬ 
cated at rtr = 10", there is an excess due to scattered light of up 
to about 0.35 (0.17) magarcsec^^ for r ^ 10". When the PSF 
is instead truncated at - 25", the corresponding convolved 
model becomes significantly fainter than the model convolved 
using the full PSF for r k, 23" (r ^ 24"). The correct amount of 
scattered light at, say, r = 40" is only achieved when the PSF 
is not truncated within, say, r - 44" (40" x 25"/23"). Com¬ 
pared to the centre, the PSF is about 12 (14) mag arcsec^^ fainter 
at r = 40". The model surface brightness is about 4.1 (2.4) 
mag arcsec^^ fainter than the convolved structure at r = 31 ”2 
(Ap ^ 14 mag arcsec^^), and the required S/N^44 (S/N^9.1). At 
r = 1975 {Ap ^ 8.7 mag arcsec^^), the corresponding values 
are about 0.78 mag arcsec^^ and S/N=s2.1 (0.26 mag arcsec^^ and 
S/N=sl.3). 

The face-on galaxy shows excess light for r k 12", where 
the break radius is at rb, = 772, Fig. A.2e. The model surface 
brightness is about 5.6 (4.1) magarcsec^^ fainter than the con¬ 
volved structure at r = 3775 {Ap ^ 14magarcsec“^); the S/N 
value that would be required to extract the intensity structure 
of the model structure is about 170 (44). At r = 26" {Ap ^ 
8.7 mag arcsec”^), the corresponding values are 1.6 mag arcsec^^ 
and S/N^4.4 (0.6mag arcsec^^ and S/Nail.7). These values are 
0.34-1.7 mag arcsec^^ brighter than in the case of the Sersic-type 
galaxy. 

For the intermediate-size edge-on galaxy model in Paper 1, the 
model surface brightness is about 6.4 (4.7) mag arcsec^^ fainter 


than the convolved stmcture at z = 2875 {Ap ^ 14 mag arcsec^^). 
At z = 19" {Ap ^ 8.7magarcsec“^), the corresponding value 
is 2.1 (0.94) magarcsec^^. Compared to the Sersic-type galaxy 
model, these values are 0.68-2.3 mag arcsec^^ brighter. 


Appendix B.3: The two targe gaiaxies 

Using either one of the two PSFs, the convolved model shows 
significant excess due to scattered light, beginning at r 160", 
Fig. A.2c. The halo is not reproduced at all when the PSF is 
truncated at - 10". The scattered-light halo of the convolved 
model is too faint for r ^ 230" when the PSF is truncated at 
r - 250". The PSFs are about 20 mag arcsec^^ fainter at r = 330" 
than at the centre. To measure the modelled intensity at r ^ 250", 
which is about 1.4 mag arcsec^^ fainter than the convolved model, 
the required S/N « 3.6, the PSF likely needs to be more accurate. 
The convolved profiles of the two PSFs nearly overlap since they 
are the same for r > 200". As in the case of the edge-on galaxy 
in Paper I, the surface brightness of the model convolved with 
PSFk 7 i in Fig. A.2c is brighter at large radii than it is in Figs. A.2a 
and A.2b; this is expected, because, compared to where r ^ 30", 
the PSF slope is shallower at larger radii. 

The face-on galaxy shows excess light for r ^ 180", where 
the break radius is at = 72", Fig. A.2f. The model surface 
brightness is about 5.0magarcsec“^ fainter than the convolved 
strucmre at r = 375" {Ap ^ 14 mag arcsec^^), the required 
S/Ns» 100. At r = 255" {Ap ^ 8.7 magarcsec^^), the correspond¬ 
ing values are 1.1 magarcsec^^ and S/N^2.8. These values are 
0.58-1.7 mag arcsec^^ brighter than in the case of the Sersic-type 
galaxy. 
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Table A.l. Observations. 


Galaxy 

Sect. 

Ref. 

Band 

Telescope / Observatory 

Date 

Exp.time 

Comments 

NGC 5907 

Paper I 

(1) 

R Harris 

Kitt Peak Nat. Obs. No. 1 0.9m 

1990-04-29, 30 

30minx6 

375, O'.'77/px 



(2) 

V 

Canada-Erance-Hawaii Tel. 

1995-06-03- 

30minx5 

< 1", O'.'206/px 





(CEHT) 

1995-06-07 





(3) 

6660 A“ 

Beijing Astr. Obs. 0.6/0.9m Tel. 

1995-01-31- 

10, 20min 

4", l'.'71/px 






1995-06-27 

xl37 

23 nights 

IC 5249 

C.l 

(4) 

Rc” 

0.6m Boiler Chivens Tel. 

1997-07-04, 

20minx7 

4" , O'.'645/px 





Mount John Univ. Obs. 

1997-07-06 



NGC 4565 

C.2 

(5) 

r 

48" Schmidt, Palomar Obs. 

1976-03-30- 

55minx2 

l",2.013"/px 






1976-04-02 




C.2 

(6) 

B, r 

48" Schmidt, Palomar Obs. 

1975-03-31- 

190min 







1977-03-13 





(7) 

V 

2.52m Nordic Optical Tel. 

1994-02-10, 11 

25min 

1", O'.'462/px 



(8) 

6660 A“ 

60/90cm Schmidt Tel., Xinlong 

1995-01-28- 

42.8h 

2'.'9^'.'4, l'.'7/px 





Stat., Nat. Astr. Obs. China 

1997-05-28 



NGC 4244 

C.3 

(5) 

r 

48" Schmidt, Palomar Obs. 

1976-03-30- 

55min 

l",2.013"/px 






1976-04-02 





(9) 

R Harris 

24/36" Burrell Schmidt Tel. 

1997-03 (5n), 

20min 

2703/px 






1997-04 (6n) 

x66 


FGC 310, 

3.2 

(10) 

R 

du Pont 2.5m, 

1997-09-21 

360s 

079, 07259/px 

FGC 1285 

3.2 

(10) 

R 

Las Campanas Obs. 

1998-04-15 

900s 

171, 07259/px 

LSBGs 

3.3 

(11) 

g, r, i 

SDSS, Apache Point 2.5m 




NGC 7280 

4.2 

(12) 

R 

2.5m Isaac Newton Tel. 

2003-09-20 

10,200s 

170-175 

NGC 4102 

4.2 

(12) 

R 

Wide-Eield Camera 

2004-03-17 

1200s 

077-374 

NGC 4477 

D.2 

(12) 

r 

SDSS, Apache Point 2.5m 




NGC 3507 

D.2 

(12) 

r 

SDSS, Apache Point 2.5m 




NGC 2880 

D.2 

(12) 

R 

Wisconsin Indiana Yale 

1997-03-02 

300s 






NOAO (WIYN) Tel. 




ESO 400-G043 

5.2.1 

(13) 

B Cousins 

ESO 2.2m Tel. 

1989 

30min 

< 1" 



(14) 

V 

ESO Multi-M. Ins. (EMMI) 

2005 

77min 






New Tech. Tel. (NTT) 




Mrk5 

5.2.2 

(15) 

B 

NOT 

1999-01 

3100s 


ESO 350-IG038 

E.l 

(13) 

B Cousins 

ESO 2.2m Tel. 

1984 

25min 

< 1" 



(16) 

y 

NOT, MOSaic CAmera 

2008 

40min 

07217/px 

Mrk 297 

E.2 

(17) 

B Mould 

Kitt Peak No. 1 0.9m Tel. 

1983-05-15-17 

50min 

171 



(15) 

B 

Calar Alto 2.2m Tel. 

1998-08 

2000s 


UM465 

E.3 

(18) 

B 

NOT, Andalucia Paint Obj. 

2004 

40min 






Spectrograph (ALPOSC) 




VCC 0001 

E.3 

(19) 

g, r, i 

SDSS, Apache Point 2.5m 




NGC 4551 

5.3 

(20) 

V 

Various telescopes 




NGC 4473 

5.4 

(20) 

V 

Various telescopes 




NGC 4374 

5.4 

(20) 

V 

Various telescopes 





Notes. Column 1, galaxy name; Col. 2, the galaxy is referred to in this section or appendix; Col. 3, reference; Col. 4, considered bandpass out of 
the used ones; Col. 5, the telescope and observatory of the measurements; Col. 6, dates of the measurements; Col. 7, total exposure time of the 
measurements; Col. 8, additional comments such as seeing and spatial resolution, where this information is available. 

The BATC filter is centred on 6660 A and has a bandwidth of 480A. ® The observations were made with a ‘step function’ filter that includes 
60 per cent of the Rc passband; the measurements were converted to the r band. The photographic exposures were on Illa-J emulsion with a 
Wratten-2C filter, which gives a bandpass of about 4000-5400A; the resulting magnitude scale is denoted J. 

References. (1) Morrison et al. (1994); (2) Lequeux et al. (1996); (3) Zheng et al. (1999); (4) ABC99; (5) van der Kruit (1979); (6) Jensen & Thuan 
(1982); (7) Naslund & Jorsater (1997); (8) W02; (9) F99; (10) DBOO; (11) BZCIO; (12) EPB08; (13) BO02; (14) M13a; (15) COl; (16) MIO; 
(17) P96; (18) M13b; (19) Meyer et al. (2014); (20) K09. 


For the large edge-on galaxy model in Paper I, the model 
snrface brightness is abont 5.8magarcsec“^ fainter than the 
convolved stmctnre at z = 285" (A/r ^ 14 mag arcsec^^). At 
z = 190" (A// ^ 8.7magarcsec“^), the corresponding valne 
is 1.7 mag arcsec^^. Compared to the Sersic-type galaxy model, 
these values are 1.2-2.5 mag arcsec^^ brighter. 


Appendix C: More on edge-on disc gaiaxies 


The three galaxies discnssed here are: IC 5249, Appendix C.l; 
NGC 4565, Appendix C.2; and NGC 4244, Appendix C.3. All 
considered existing and derived model parameters of these galax¬ 
ies as well as the two thick-disc galaxies FGC 310 and FGC 1285 
(Sect. 3.2) are given in Table C.l. 
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Fig. A.2. Surface-brightness profiles of example models of three Sersic-type galaxies (the upper panels a-c) and three face-on disc galaxies (the 
lower panels d-f) that illustrate effects of a radial truncation of PSFkti and PSFv,om- The three columns show: a) and d) two small galaxies, b) and 
e) two intermediate-size galaxies, and c) and f) two large galaxies. In each panel, the galaxy model profile is drawn with a thick black solid line, and 
the model profiles that were convolved with PSFkti (PSFv/om) are drawn with thin purple (blue) lines. The profile of the model that was convolved 
with the complete PSF is drawn with a thin solid line. Profiles are also shown where each PSF was truncated at rtr.i iraa) with a dotted (dash-dotted) 
line in each panel; a) and d) i = 4" and r „.2 = 10"; b) and e) = 10" and ri ^2 = 25"; c) and f) r,rj = 10" (lines fall atop the model lines and 
are not visible) and = 250" (Table A.2). The dash-tripple-dotted lines show the result of one-dimensional models that were convolved with the 
respective one-dimensional PSF. 


Appendix C. 1: The halo of 1C 5249 

The outcome of f?-band observations of the edge-on Sc/Sd-type 
galaxy IC 5249 are presented by ABC99, who measure PSFa 99 
that extends to r = 39". They derive the PSF following the 
approach of Morrison et al. (1994); this /?-band PSF is brighter 
than other typical PSFs, see Fig. A.l. In the outer region, the PSF 
agrees reasonably well with both PSF,_ 3 n, and PSF,_om- ABC99 
study the contribution of the PSF to the halo by convolving it with 
an image of the galaxy truncated at z = 5'.'7. The resulting surface- 
brightness profiles are significantly fainter than the measurements 
(see their fig. 6 ), and they conclude that the halo is real. 

The authors model the surface-brightness profiles of IC 5249 
by superposing a thick disc on a thin disc. I modelled IC 5249 
with one single thin disc with similar, but not identical parameters 
as ABC99, where zo — 4'.'7, - 63", and - 1.1 convolved 

the /?-band model with the usual three PSFs that I use in this 
paper, but also with PSFa 99 and PSFj^m. Before the convolution, 
I extrapolated PSFa 99 with an power-law profile where r > 
40", cf. Fig. A.l. Surface-brightness profiles of both models and 
measurements are shown in Fig. C.la, along vertical-axis profiles 


of the part of the galaxy that ABC99 refer to as region B, which 
is the region most closely centred on the minor axis. 

Because of the larger distance to IC 5249, the surface- 
brightness profiles are steeper than in NGC 5907 (where zo - 
15"). Intensity and colour profiles are thereby also more sensitive 
to varying PSFs in IC 5249. With PSFyom, the scattered-light 
halo begins already at z 8 ". The /?-band model convolved with 
PSFa 99 agrees slightly better with the measurements than the 
model convolved with PSF, 3 ni; compare the two measurements at 
r ^ 17" and r ^ 20". PSFa 99 extends to 3876, which is enough 
to analyse measurements out to about r - 25" (3876/1.5). The 
six measurements for r ^ 22 " are fainter than either model, they 
are too low for either PSF and have too small upper error limits. 
The halo is about 0.5 mag arcsec^^ too bright with PSF,- 3 ^- No¬ 
tably, the required accuracy at the location of the three outermost 
measurements, which would allow the data to be deconvolved to 
an exponential profile (ignoring the achievable accuracy of the 
sky subtraction), are about S/N=40, 300, and 500, respectively. 

The colour profiles show red excess for z > rno = 8 ", 
Fig. C.lb. The excess at z = 20" is A{R-i) — 0.85- 
1.35 mag arcsec^^, compared to the modelled profile, which 
grows to A(R-i) - 1.35-2.00magarcsec^^ at z = 30". The 
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Table C.l. Model parameters of the 7 edge-on disc galaxies considered and referred to here. The distances are in this work only used to scale the 
upper x-axis in the surface-brightness profile plots. Details of the original observations are given in Table A.l. 


Galaxy 

D 

Band 

ris 

single/thin disc 


thick disc 


bulge 


Ref. 

Fig. 

hr 

Zo 

Bo,o 

hr 

Zo 

bq,q 

n 

re 

Ap 

NGC 5907 

11 

R Harris 

2 

90 

16 

- 

- 

- 

- 

- 

- 

- 

1 



11.7 


2 

90 

15 

- 

- 

- 

- 

- 

- 

- 

2 


IC 5249 

36 

Rc 

CO 

63 

3.3 

- 

34 

9.7 

- 

- 

- 

- 

3 

C.l 


36 


1 

63 

4.7 

- 

- 

- 

- 

- 

- 

- 


C.l 

NGC 4565 

10 

J 

1 

110 

16 

20.1 

- 

- 

- 

- 

- 

- 

4,5 

C.2 


10 

B 

1 

124 

17 

- 

120 

24 

- 

- 

- 

- 

6, T 

C.2 



r 

1 

100 

17 

- 

110 

26 

- 

- 

- 

- 

6, 7" 

C.2 


10 

V 

1 

120 

12 

20.2 

120 

27 

21.8 

- 

- 

- 

8^ 

C.2 


14.5 

6660A 

1 

115 

17 

22.32 

160 

36 

25.52 

- 

- 

- 

9 

- 


10 


1 

113 

20.6 

- 

- 

- 

- 

1 

3.0 

-0.7 


C.2 

NGC 4244 

5 

J 

1 

no 

24 

20.6 

- 

- 

- 

- 

- 

- 

4,5 

C.3 


3.6 

R 

oo 

105.4 

14.1 

20.47 

- 

- 

- 

- 

- 

- 

10 

C.3 


3.6 


2 

105 

13 

- 

- 

- 

- 

- 

- 

- 


C.3 

FGC310" 

80.6 

R 

1 

8.7 

1.95 

21.19 

- 

- 

- 

- 

- 

- 

11 

2a&2b 


80.6 

R 

1 

8.4 

1.6 

21.3 

9.6 

2.9 

22.5 

- 

- 

- 

11 

2a&2b 


80.6 


1 

8.4 

1.6 

21.3 

- 

- 

- 

- 

- 

- 


2a&2b 

FGC 1285" 

18.8 

R 

1 

19.7 

6.63 

20.99 

- 

- 

- 

- 

- 

- 

11 

2e&2f 


18.8 

R 

1 

17.3 

4.4 

21.2 

23.6 

10.1 

22.2 

- 

- 

- 

11 

2e&2f 


18.8 


1 

17.3 

9.6 

19.5 

- 

- 

- 

1 

4.4 

-1.6 


2e&2f 


Notes. Column 1, galaxy name; Col. 2, distance (Mpc); Col. 3, bandpass used in the observations; Col. 4, the «s-parameter in Eq. (4); Cols. 5-7, 
single-disc or thin-disc scale length (arcsec), scale height (arcsec), and intensity at the object centre (mag arcsec^^); Cols. 8-10, thick-disc scale 
length (arcsec), scale height (arcsec), and intensity at the object centre (mag arcsec^^); Cols. 11-13, bulge Sersic index, half-light radius (arcsec), 
and magnitude difference A/r = /jq.o - Col. 14, references, rows without a number refer to this study; and Col. 15, results are shown in this figure. 

The values are taken from the sech^-sech^-exp-model in table 6 of (9). ® The authors fit the galaxy profiles with three discs where the parameters 
of the third disc are = 200" and zo = 7". Parameter errors are omitted, considering the uncertainty due to scattered light that appears to affect 
the models of (13). 

References. (1) Morrison et al. (1994); (2) Paper I; (3) ABC99; (4) Jensen & Thuan (1982); (5) Shaw & Gilmore (1989); (6) Naslund & Jorsater 
(1997); (7) W02; (8) F99; (9) YD06. 


smaller (larger) values used PSFj/ 3 n, and PSF, 3 n) (PSFy,om and 
PSF,;om)- As is the case with the separate surface-brightness pro¬ 
files, the large range of magnitudes between the two colour pro¬ 
files indicates a strong sensitivity to the PSF. 

IC 5249 is a good example where it is clear that the observed 
halo is strongly affected by scattered light. It is possible to draw 
this conclusion because of the availability of a PSF with unusually 
bright wings observed together with the object. NGC 3957 is 
another galaxy with surface-brightness structures that are fitted 
with a somewhat similar scale height as IC 5249 (Jablonka et al. 
2010). Whilst similar models also appear to explain its halo as 
scattered light (not shown here), a bar-like region at intermediate 
radii (Pohlen et al. 2004) complicate the analysis. 

Appendix C.2: The thick disc and the hate of NGC 4565 

The galaxy NGC 4565 is (along with NGC 891) an outstanding 
example of a bright edge-on galaxy similar to the Milky Way 
(van der Kruit 1984). van der Kruit (1979) observed NGC 4565 
in the J band, and the resulting plates were analysed by KS81a, 
who smooth the data and find a thick disc. Later, Jensen & Thuan 
(1982) present profiles in the B band and the r band of new 
observations (their table 10 and fig. 18), which they compare with 
previous measurements of the same galaxy (their section Ille). 
The two studies present surface-brightness profiles that differ in 


the centre parts, this is likely a consequence of how they smooth 
their data. Additional profiles are presented by Naslund & Jorsater 
(1997, V band) and (Wu et al. 2002, hereafter W02) (6660A 
band), where the latter authors find both a thick disc and a halo in 
their deep observations. W02, furthermore, claim they measure 
a PSF out to r = 1700", but it is not shown, and it is used to 
convolve a model instead of deconvolving it. The PSF is said to 
be derived from separate observations of bright stars, whilst the 
galaxy observations were made in a time span of more than two 
years; there is no information on when those observations were 
made, and the PSF data are meanwhile lost (Wu, priv. comm. 
2014). 


A study of near-infrared (NIR) observations of NGC 4565 
(3.5-5;um) from outside of the Earth atmosphere finds no proof 
for a halo (Uemizu et al. 1998). Unlike other studies, these authors 
measure a PSF with a full-width-at-half-maximum of about 25" 
out to r = 10', where it is 12.5 magarcsec^^ fainter than at the 
centre (the PSF is not shown). They use the PSF to show that 
their halo is scattered light. A quick check in fig. 1 in Paper I 
shows that the corresponding value of optical wavelength-range 
PSFs at the same radius is some 18 mag arcsec^^ fainter than at 
the centre; because of the steeper decline of the PSF in the optical 
wavelength range, one might be led to conclude that scattered 
light is of little concern in this case. 
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Fig. C.l. Minor-axis i?-band and i-band surface- 
brightness profiles versus the vertical distance 
z of models and measurements of the edge-on 
galaxy IC 5249. Panels, lines, and coloured 
bullets are as described in Fig. 2. The f?-band 
model was additionally convolved with PSFaqs 
(PSF, 3 n,) to produce the cyan-coloured solid 
(dash-dotted) line with a black border in panel 
a). Crosses, circles, and error bars in the same 
panel show region-B measurements for cuts on 
either side of the disc of ABC99. 


KS8 la fit the brightest parts of the spatially smoothed surface- 
brightness structure. Shaw & Gilmore (1989) make multiple two- 
and three-component fits of the surface-brightness structure (us¬ 
ing the data of Jensen & Thuan 1982). The final r-band and 
B-band models (their table 6) contain a thin disc and a thick disc, 
and a halo with an exponential fall-off {h - 53")- Naslund & 
Jorsater (1997) fit their V-band galaxy profile with three discs. 
W02 fit a thin disc, a thick disc and a power-law halo to their 
surface-brightness profile. They compare their fit with nine previ¬ 
ous studies. 

I fitted the measurements of Jensen & Thuan (1982) with a 
small bulge superposed on a single thin disc, where - 113", 
Zo = 2076, tis - I, n — I, re - 3", and jUoo - Be = 
-0.7 magarcsec^^. The resulting minor-axis surface-brightness 
profiles and the measurements of KS81a, Jensen & Thuan (1982), 
and Naslund & Jorsater (1997) are shown in Fig. C.2a (the mea¬ 
surements of W02 are left out as they would overfill the plot). 
The inclusion of the bulge has minor effects on the outer parts of 
the convolved profiles. 

Despite that the apparent size of NGC 4565 is much larger 
than IC 5249 and NGC 3957, the colour profiles are similarly 
affected by scattered light. Fig. C.2b. The red excess amounts 
to A (B - i) ^ 0.85-0.96 mag arcsec^^ at z = 80" and about 1.2- 
1.3 mag arcsec^^ at z = 120", where the larger (smaller) values 
were calculated using PSFyon, and PSF;_om (PSFyjm and PSF,_ 3 n,)- 


The presented profiles of NGC 4565 provide strong hints that 
scattered light can induce a structure that appears as both a thick 
disc and a halo around this galaxy. The four measured profiles are 
suspiciously similar (as is the fifth one of W02); in comparison, 
the measured profiles of NGC 5907 show larger variation. The 
studies do not provide enough information to exclude that they 
have matched their sky subtraction to achieve the very similar 
profiles of NGC 4565. To correctly deconvolve observations of 
this galaxy out to r = 100", it is necessary to use a PSF that 
extends out to, at least, r = 150" in each considered bandpass. In 
this context, it is worth to mention that the PSF is poorly known 
at large radii, and at the time of the observations it could very 
well be brighter than the lower-limit PSFkvi, and even as bright 
as PSFy 3 ni for r ^ 100". At r = 100", the assumed exponential 
model structure is about 2 mag arcsec^^ fainter than the observed 
structure. If the data would be deconvolved to recover the expo¬ 
nential profile, the required signal-to-noise is about S/N= 6.5 at 
the same radius. 

Appendix C.3: The existence of a thick disc and a haio in 
visuai waveiength-range observations of NGC 4244 

Using the photographic observations of NGC 4244 made by van 
der Kruit (1979), KS81a analyse and smooth the data, but do not 
find a thick disc. Surface-brightness profiles of more recent obser¬ 
vations in the R band are presented by F99. The authors follow 
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Fig. C.2. Minor-axis surface-brightness profiles 
versus the vertical distance z for the edge-on disc 
galaxy NGC 4565. Panels, lines, and coloured 
bullets are as described in Fig. 2. The bullets 
and circles (diamonds; orange stars) in panel a 
show measured r-band and S-band (T-band; V- 
band) values of Jensen & Thuan (1982, KS81a; 
Naslund & Jdrsater 1997). 



z (arcsec) 


the approach of Morrison et al. (1994) and measure PSFpgg that 
extends out to r = 305"; no information is provided regarding the 
saturated stars used to derive the outer profile for r S 20". PSFpgg 
is shown in Fig. A. 1 where it is seen to be similar to PSFmbh and 
PSF^gjj. Any influence of scattered light is dismissed; the PSF 
is merely scaled to the measurements of the centre intensity of 
NGC 4244 to show that the PSF is fainter at large radii. 

F99 fit a single exponential disc to the data. I also fitted the 
measurements with a single disc, where - 105", zo - 13", 
and Hs = 2.1 set the model intensity constant for z < 29". Minor- 
axis surface-brightness profiles and the measurements of both 
KS81a and F99 are shown in Fig. C.3a; only five error bars were 
extracted from fig. 9 in the latter reference since the figure is so 
small (some of the data points are also reproduced in a larger 
figure in Morrison 1999). Tikhonov (fe Galazutdinova (2005) and 
Seth et al. (2005) use HST data and count resolved red-giant 
stars out to large vertical distances. The re-analysed and scaled 
data points of the second data set (Streich et al. in prep.) are 
shown together with the photometry data in Fig. C.3a. The data 
of Tikhonov & Galazutdinova (2005) were scaled to match the 
second data set where r ^ 80". 

The colour profiles in Fig. C.3b show red excess for z ^ 45", 
which amounts to A (/? - i) ^ 0.36-0.39 mag arcsec^^ at z = 80" 
and about 0.75-0.83 mag arcsec^^ at z = 120". This galaxy is 
intermediate in apparent size between NGC 3957 and NGC 4565, 


and despite a reported lack of a thick disc or halo, the models 
reveal significant amounts of red excess. 

Both sets of photometry measurements of NGC 4244 show 
good agreement with the surface-brightness profiles of the con¬ 
volved models. As with the other edge-on galaxies, photometry 
of this object also reveals a scattered-light halo, which deviates 
from an exponential profile for r k, 80". This galaxy is somewhat 
fainter than NGC 5907, which could explain why there are fewer 
measurements of the scattered-light halo at large radii. The star- 
count profiles deviate from the photometry profiles, possibly due 
to how stellar populations are calculated from the red-giant star 
counts. The important attribute is that the single-star counts result 
in a halo that appears to be some 2 mag arcsec^^ fainter than the 
photometry measurements for r k 100". 


Appendix D: More on face-on disc gaiaxies 

Appendix D. 1: Model parameters 

The model parameters of all six face-on disc galaxies considered 
and referred to here and in Sect. 4 are given in Table D. 1. 
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Fig. C.3. Minor-axis surface-brightness pro¬ 
files versus the vertical distance z for the edge- 
on disc galaxy NGC 4244. Panels, lines, and 
coloured bullets are as described in Fig. 2. The 
black bullets and error bars (diamonds) in panel 
a show measured /?-band (T-band) values of F99 
(KS81a). The grey boxes and stars show scaled 
single-star measurements derived from HST star 
counts of Tikhonov & Galazutdinova (2005) and 
Streich et al. (in prep.), respectively; the val¬ 
ues of the latter reference are only shown for 
z > 40". The best-fit profile of F99 (zo = 14'.T) 
is shown with a black dotted line. 


Table D.l. Fitted f?-band model parameters of the 7 face-on disc galaxies considered and referred to here. The galaxy type, distance, inclination, 
and the disc parameters are identical to the values of EPB08 (their tables 1 and 4 and fig. 14). A bulge was also added to each model. Details of the 
original observations are given in Table A. 1. 


Galaxy 

Type 

D 

1 

Inner and outer discs 

hi ho yUo.i 

/^0,o 

fbi- 

Bulge 

n re 

Pc 

^abr 

Fig. 

NGC 7280 

II.o-OLR(?) [-F 111?] 

24.3 

50 

25.4 

11.2 

20.20 

17.90 

42.5 

4.0 

2.5 

17 

100 

4a, b 

NGC 4102 

II.o-OLR -H III-s(?) 

14.4 

55 

344.1 

15.4 

20.58 

17.59 

42.9 

1.0 

7.5 

17.6 

100 

4c, d 

NGC 4477 

I [or III?] 

16.5 

33 

35.7 

- 

19.87 

- 

- 

4.0 

7.0 

17.5 

200 

D.la,b 

NGC 3507 

II.o-OLR(?) [-F III?] 

14.2 

27 

36.7 

18.9 

20.31 

18.22 

75.1 

4.0 

1.5 

17.3 

130 

D.lc, d 

NGC 2880'' 

III-s 

21.9 

52 

20.0 

- 

19.77 

- 

- 

2.0 

3.0 

15.85 

120 

D.le, f 


Notes. Column 1, galaxy name; Col. 2, galaxy type; Col. 3, distance (Mpc); Col. 4, disc inclination (°); Cols. 5-9, disc parameters of EPB08, 
the scale lengths of the inner and outer discs (arcsec), the inner and outer disc centre intensities (mag arcsec^^), and the break radius between the 
inner and the outer discs (arcsec); Cols. 10-12, bulge Sersic-profile index, half-light radius (arcsec), and centre intensity (magarcsec^-); Col. 13, 
approximate anti-truncation break radius (arcsec); and Col. 14, galaxy model is shown in these figures. 

This value results in a large offset between the inner and the outer discs; I used the value po o = 17.27 magarcsec^^ instead. 


Appendix D.2: Three examples where the PSF may explain 
faint structures: NGC 4477, NGC 3507, and NGC 2880 

Surface-brightness profiles of models of NGC 4477 are shown 
in Fig. D.l a, together with the SDSS r-band measurements of 
EPB08. The profile is classified as type I or possible type 111. An 
anti-truncation break appears at r^br - 200". The scattered-light 
halo radii are rno.PSPK,! = '‘iio.PSFv.on ~ 170" and riio,PSF,.o™ ~ 


100". This galaxy is so large that the convolved /?-band profiles 
overlap as the outermost parts of the PSFs all share the same 
values. The measurements are brighter than the convolved model 
profiles for radii r ^ 210", which is only slightly larger than 
the anti-truncation radius rabr - 200". The colour profiles in 
Fig. D.lb show smaller amounts of red excess than in the other 
face-on galaxies of this study, because of the large scale length; 
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Fig. D.l. Radial R-band and i-band surface-brightness profiles versus radius for three face-on disc galaxies: a) and b) NGC 4477, c) and d) NGC 
3507, and e) and f) NGC 2880. The panels, lines, coloured bullets, and offsets are as described in Figs. 2 and 4. Diamonds (crosses) show R-band 
SDSS (INT) measurements (EPB08, their fig. 14); red symbols indicate values that are fainter than their sky-uncertainty limit. 


shortwards of r ^ 25" there is also some blue excess. The values 
are A (R - i) ^ 0.086-0.095 mag arcsec^^ at r = 100" and about 
0.10-0.24 mag arcsec^^ at r^br- 

Model surface-brightness profiles of NGC 3507 are shown in 
Fig. D.lc, also in this case with SDSS r-band measurements of 
EPB08. The anti-truncation radius =» 135" agrees with the 
scattered-light halo radius rno - 125-135". This galaxy is not 
as extended as NGC 4477, but the convolved profiles that were 
calculated using PSFkvi and PSFyom overlap for r > rno- The 
measurements are brighter than the profiles for r > rabr - 130" 
also here. The colour profiles in Fig. D.ld are reminicent of 
those for NGC 4102 in Fig. 4d. There is red excess at all radii 
outside the break radius rbr = 75'.'1, as with the other galaxies 
with such a break. At larger radii, the colour excess values are 
A (R - i) ^ 0.11-0.13 mag arcsec^^ at r = 100" and about 0.25- 
0.30 mag arcsec^^ at rabr- The two radii rno and rabi are very 
similar for both NGC 4477 and NGC 3507, and it cannot be 
excluded that more extended and brighter SDSS PSFs could 
explain the suggested type III profile in the outermost parts of 
these two galaxies. 

The final galaxy considered is NGC 2880, which surface- 
brightness profiles are shown in Fig. D.le, together with the 
measurements of the WIYN telescope. The parameters of this 
galaxy are similar to those of NGC 3507, but no inner break radius 
is used with this galaxy. The scattered-light halo radii ri m are 
again similar to the anti-tmncation radius rabr = 120". The PSFs, 
and PSFv, 3 in in particular, seem to be able to explain the measured 
profile out to r = 150". There are no known measurements of the 
PSF of the WIYN telescope at larger radii. The colour profiles 
in Fig. D.lf show red excess outside of r ^ 25". The values 
are A(R-i) ^ 0.13-0.17 mag arcsec^^ at r = 50" and about 
0.34-0.36 mag arcsec^^ at rabr- 


Appendix E: More on BCGs 

Three BCGs are studied in the following subsections: ESO 350- 
IG038, Appendix E.l; Mrk 297, Appendix E.2; and UM 465, 
Appendix E.3- The model parameters of all nine BCGs considered 
and referred to here and in Sect. 5.2 are given in Table E.l. 

Appendix E. 1: The extreme BCG ESO 350-IG038 (Haro 11) 

Bd02 present surface-brightness profiles of observations in the B 
band. The observed surface-brightness profiles were fitted using 
ellipses with the ellipticity e = 0.20 (as Micheva et al. 2010 
report later, hereafter Ml0). The latter authors present more recent 
profiles of observations in the V band. In this case, the measured 
ellipticity is e = 0.20-0.32. 

The first authors fit the surface-brightness profile for a part 
of the radial range with a Sersic profile. They also fit the mea¬ 
surements in the same radial range with an exponential disc. I 
matched the observed R-band measurements with a convolved 
model where I used two Sersic profiles: a bulge where n = 1.0, 
Te = 077, and = 13.5 mag arcsec^^, and a host where 
n - 1.0, Te = 570, and - 22.0 mag arcsec^^. Additionally, 
Re.B - Re,I = 0.5 mag arcsec^^. To match the measurements, the 
profile is truncated for r 276. The surface-brightness profiles 
and the two sets of measurements are shown in Eig. E.l a. It is dif¬ 
ficult to fit the measurements well at all radii, both the innermost 
and the outermost measurements are poorly fitted in either case. 
The convolved profiles are up to 1.0 mag arcsec^^ too bright for 
r < 270. The measurements appear to be fairly well matched with 
a real PSE that is similar to PSEyom- The more recent V-band 
measurements were fitted with ellipses of different ellipticity, and 
appear radially offset at the centre of the object - they are poorly 
fitted in any case. 
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Table E.l. Model parameters of the 9 BCGs that are considered here. Only those values are shown where a value is provided. The distances are in 
this work only used to scale the upper x-axis in the surface-brightness profile plots. Details of the original observations are given in Table A.l. 


Galaxy 

D 

Band 

Low surface-brightness component 

Starburst component 

rx 

Ref. 

Fig. 

n 

fe 

h 

B 

n 

fe 

h 

Ap 

ESO 400-G043 

77 

B 

17.1 

1.13 mas 

- 

-38.85 

- 

- 

- 

- 

- 

1“ 

6 


77 

B 

1 

11 

6.7 

23.1 

- 

- 

- 

- 

- 

1“ 

6 


79.2 

V 

1 

7.21 

4.31 

20.64 

- 

- 

- 

- 

- 

2* 

6 


79.2 

V 

1 

15.7 

9.39 

23.60 

- 

- 

- 

- 

- 

2“^ 

6 


79.2 


1 

4.3 

2.6 

- 

1 

1.5 

0.90 

4.75 

5.0 


6 

Mrk5 

13.3 

B 

1 

8.02 

4.81 

21.09 

- 

- 

- 

- 

- 

3'' 

8 

UGCA 130 

13.3 

B 

2.83 

3.35 

- 

20.23 

- 

- 

- 

- 

- 

4 

8 


13.96 

B 

0.99 

9.6 

- 

21.18 

- 

- 

- 

- 

- 

5 

8 


13.3 


1 

8.0 

4.8 

- 

1 

2.0 

1.2 

0.7 

- 


8 

ESO 350-IG038 

80 

B 

19.8 

0.080 mas 

- 

-53.95 

- 

- 

- 

- 

- 

r 

E.l 

Haro 11 

80 

B 

1 

14 

8.3 

23.5 

- 

- 

- 

- 

- 

r 

E.l 


82 


1 

5.0 

3.0 

- 

1 

0.7 

0.42 

8.50 

2.5 


E.l 


82 


11.0 

8.25 mas 

- 

-41.20 

- 

- 

- 

- 

- 


E.4 

Mrk 2977 

63 

B 

1 

16.3 

9.72 

22.43 

- 

- 

- 

- 

- 

6 

E.2 

NGC 6052/54 

64.3 

B 

1 

13.3 

7.96 

22.30 

- 

- 

- 

- 

- 

3« 

E.2 

UGC 10182 

65.10 

B 

0.84 

9.25 

- 

21.35 

- 

- 

- 

- 

- 

5 

E.2 


64.3 


5.5,: 

3.5, 18.48, 16.6, 

16, 22* 


1 

2.0 

1.2 

- 

- 


E.2 

UM 465 

20.7 

B 

1 

10.8 

6.48 

19.97 

- 

- 

- 

- 

- 

7* 

E.3 

UGC 6877 

20.7 

B 

1 

15.7 

9.37 

21.82 

- 

- 

- 

- 

- 

7“^ 

E.3 

Mrk 1308 

20.7 


1 

9.53 

5.70 

- 

1 

2.5 

1.5 

2.4 

- 


E.3 

VCC 0001 

32 

r 

1 

2.6 

1.5 

20.3 

- 

- 

- 

- 

- 

8 

E.5 


32 


1 

4.9 

2.9 

- 

- 

- 

- 

- 

- 


E.5 


Notes. Column 1, galaxy name, some alternative names are provided indented on subsequent lines; Col. 2, distance (Mpc); Col. 3, bandpass used in 
the observations. Cols. 4-7, Sersic parameters of the low surface-brightness component in the form of the index, the half-light radius (arcsec), the 
corresponding scale length when « = 1 (arcsec), and the intensity at r = 0 (n = 1) or (mag arcsec^^, n t 1); Cols. 8-11, Sersic parameters of the 
starburst component in the form of the index, the half-light radius (arcsec), the corresponding scale length when n = 1 (arcsec), and the magnitude 
difference Ap = p - (mag arcsec^^); Col. 12, inner truncation radius (arcsec); Col. 13, references, rows without a number refer to this study; and 
Col. 14, results are shown in this figure. 

The fit is based on measurements in the radial range 16 < r < 32". ® The scale length and the central surface brightness are based on 
measurements in the radial range where 24 < /r < 26 mag arcsec^^. The scale length and the central surface brightness are based on measurements 
in the radial range where 26 < yu < 28 mag arcsec^^. The fit is based on measurements in the radial range 13 < r < 32". The fit is based on 
measurements in the radial range 20'.'6 < r < 41'.'3 (8-16 kpc). Besides the low surface-brightness component, the parameters of the two inner 
components are given in Sect. E.2. The fit is based on measurements in the radial range 30'.'3 < r < 36'.'5. *** The six values specify the five 
face-on galaxy parameters hi (arcsec), ha (arcsec), /ro.i (mag arcsec^^), yuo.o (mag arcsec^^), and rbr (arcsec), as well as p^ of the bulge (mag arcsec^^). 

References. (1) BO02; (2) M13a; (3) COl; (4) C05; (5) Amortn et al. (2009); (6) P96; (7) M13b; (8) Meyer et al. (2014). 


The three B - I colour profiles are shown in Fig. E.lb. The 
red excess at r = 20" is A(B - /) ^ 0.70 mag arcsec^^ with 
PSFv;om and PSF, om, and about 0.35 mag arcsec^^ with PSFy^m 
and PSF, 3 ni. The corresponding values at r = 35" are A(B - /) ^ 
1.7 mag arcsec^^ and about 1.0 mag arcsec^^, respectively, which 
are very similar to the corresponding values for ESO 400-G043 
at r = 30". 

As in the case of ESO 400-G043, which was modelled using 
similar small-valued parameters, the scattered-light halo of ESO 
350-1G038 begins already at small radii. The amount of red excess 
at larger radii is therefore determined by the shape of the B-band 
and the /-band PSEs. As with the other studied galaxies in this 
paper, a hypothetically temporally varying PSF can itself explain 
why the more recent V-band observations of MIO simultaneously 
show more scattered light and less red excess in the halo (they 
use the colour V - K). 

There are several separated bright starburst regions at the cen¬ 
tre of this BCG; their placement suggests the inapplicability of 


a symmetric one-dimensional function to describe the co-added 
structure of scattered light. The measurements of MIO are brighter 
than the B-band measurements, suggesting possible differences 
between how the asymmetric surface-brightness structure is sam¬ 
pled, and the PSF could also be different. The goal here was 
not to reproduce the true model profile, but to illustrate that the 
geometrical properties of this galaxy suggest that scattered-light 
effects dominate the measured halo. 

To demonstrate the importance of an accurately determined 
PSF, I calculated a set of alternative models where I set n = 11.0, 
Te = 8.2 mas, and - -41.2 mag arcsec^^. The profiles are again 
truncated for r $ 276. The resulting surface-brightness profiles 
are shown in (the online-only) Fig. E.4a. The convolved profiles 
of these alternative models and the models in Fig. E. la agree fairly 
well with the B-band measurements. There is, in this case, slightly 
less red excess than before, compare Eigs. E.4b and E.lb. Here, 
the red excess at r = 30" is A(B - /) = 1.0 mag arcsec^^ with 
PSEv,om and PSE, om, and about 0.60 mag arcsec^^ with PSEy^m 
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Fig.E.l. Surface-brightness profiles versus ra¬ 
dius for the BCG ESO 350-IG038 (Haro 11). 
Panels, lines, and coloured bullets are as de¬ 
scribed in Fig. 2. The /-band profiles are offset 
by -1.5 mag arcsec^^ from the //-band profiles 
for increased visibility. Grey and black bullets 
and error bars mark measured B-band (BO02) 
and T-band (MIO) values, respectively. The grey 
dotted (dashed) line shows the best-fit Sersic- 
type model (exponential fit to the disc in the 
range 20 < r < 41") of BO02. 


and PSF,_ 3 n,. Since the exact shape of the PSF is unknown and 
because of the asymmetry due to the brightest (starburst) regions 
of the object, it is difficult to deconvolve the provided surface- 
brightness structure to get the real physical structure. 

There seem to be three possible scenarios regarding the host 
of ESO 350-IG038: there is no host (unlikely, but not excluded 
with the current data), the host is of disc type, and the host is more 
massive (n ^ 10). It is impossible to say which one is correct 
before the scattered light of the centre regions is first removed. 
There is no strong physical basis for either of the two profiles 
that BO02 measure in the outer region. The galaxy should be 
analysed after the observations are deconvolved - accounting for 
the maximum possible useful radius when considering the sky 
subtraction and the PSF. 


Appendix E.2: The BCG Mrk 297 

P96 present B-band profiles of Mrk 297. The authors model the 
surface-brightness profile with three components: a low surface- 
brightness region Ie, a plateau Ip, and an inner part /q. These 


components are defined as follows 
h = 7e,o exp 

Ig = 7g,o exp 


r 

aJ 


-0.51 - 
7 


, h = 7p,o exp 
2 


and 


(E.l) 

(E.2) 


where the surface-brightness profile of Mrk 297 is fitted with a - 
9772, = 1370, Tj = 2.46, j = 3737, /e,o = 22.43 mag arcsec-^, 

Ipo = 20.67 magarcsec^^, and /g,o = 20.21 magarcsec^^. With 
the exception of the low surface-brightness component Ip, the 
values of this fit are omitted in Table E. 1. COl present more recent 
observations in the B band and also fit an exponential profile in a 
small part of the radial extent. Amorm et al. (2009) use the same 
data and make a two-dimensional Sersic-profile fit of the stellar 
host. 

I defined a bulge and face-on disc galaxy model to get a 
reasonable agreement between convolved structures and the mea¬ 
surements of P96 and COl. The parameters of the bulge were 
n - - 270, and pe = 22 mag arcsec^^, and those of the face- 

on galaxy were hi - 575, ho — 375, po,i - 18.48 mag arcsec^^. 
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Fig. E.2. Surface-brightness profiles versus ra¬ 
dius for the BCG Mrk297. Panels, lines, and 
coloured bullets are as described in Fig. 2. The 
/-band profiles are offset by -1.5 magarcsec^^ 
from the 6-band profiles for increased visibil¬ 
ity. Grey (black) diamonds and error bars mark 
measured 6-band values of COl (P96). The grey 
(black) dotted line shows the best-fit Sersic-type 
model of COl (P96). The dash-tripple-dotted line 
shows the fitted profile of Amorin et al. (2009). 


/*o,o = 16.6 mag arcsec^^, and - 16". I also set B - fiej - 
0.5 mag arcsec^^. Surface-brightness profiles of the models and 
the three sets of profile fits are shown in Fig. E.2a. 

The scattered-light halo radius rno - 22" for PSFyom and 
about 30" for PSFkti- The measurements of both P96 and COl 
overlap the models well in all radial regions; the values deviate 
slightly for r ^ 5", but that deviation is unimportant to the discus¬ 
sion here. Each study calculated a fit to the host component using 
the outermost measurements. Their scale lengths differ slightly 
as the outermost measurements that are used in the fit differ. The 
profiles of this smdy suggest that the difference in the outermost 
region occurs as the PSE of the two separate observations were 
different, and both fits merely measure the scattered-light struc¬ 
ture for r ^ 20". Another possibility is that the accuracy of the 
sky subtraction differs (see below). The PSE at the time of the ob¬ 
servations seems to be closer to PSEyom with the measurements 
of COl, whilst that of the measurements of P96 was slightly 
brighter, but not as bright as PSEy^m. The profile is reminiscent 
of a type 111-s structure as is seen in face-on disc galaxies. The fit 
of Amorfn et al. (2009) is more difficult to interpret as it differs 
from both sets of measurements outwards of r ^ 20 ". 

Despite the small value on B - fiej, the B - I colour 
profiles in Eig. E.2b show dramatic amounts of red excess 
for r ^ 18". Eor example, at r = 30" the red excess is 
A(B - I) ^ 0.75 magarcsec^^ with PSEyom and PSE,_om, and 


about 0.59 mag arcsec^^ with PSFy^m and PSF, 3 ^- The values 
increase dramatically at larger radii. 

It cannot be excluded that the measurements of P96 and COl 
differ because of differences in the sky subtraction. The same 
result is achieved with models of Mrk 35 and Mrk 314 (not 
shown here). The measurements are ambiguous when both sky 
subtraction and scattered light are considered. 

Appendix E.3: The BCG UM 465 

M13b present deep surface-brightness profiles for 21 BCGs. Sev¬ 
eral of the BCGs are too small, too large, show little intensity 
excess, or are severely asymmetric. These are all factors that com¬ 
plicate an analysis of the influence of scattered light as presented 
here. I considered the following five BCGs to be more easy to 
analyse with the high-seeing PSEs I use, as they neither seem to 
require a deconvolution of the images nor extend to large radii 
were the intensity of the PSE is unknown: UM 465, UM 483, UM 
491, UM 538, and UM 559. Out of these BCGs, I chose to scruti¬ 
nize the surface-brightness profiles of the most symmetric galaxy 
of these, UM 465. The authors point out that the 6 -band PSE for 
this galaxy is different from the other filters in the centre; there 
is no other information on this issue. The ellipticity is measured 
as e = 0.16. The authors fit their measured surface-brightness 
profiles with two exponential discs in two separate magnitude 
ranges. 
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Fig. E.3. Surface-brightness profiles versus ra¬ 
dius for the BCG UM 465. Panels, lines, and 
coloured bullets are as described in Fig. 2. The 
/-band profiles are offset by -1.5 magarcsec^^ 
from the S-band profiles for increased visibility. 
Diamonds and error bars mark measured //-band 
values of Ml3b. The dotted (dashed) line shows 
the best-fit exponential model fit of Ml3b for 
24 < fiB < 26 (26 < < 28) mag arcsec^^. 


I chose the model parameters of a bulge where n = 1.0 
and Te = 2"5, which is superposed on an exponential-disc host 
where n - 1.0, - 9'.'53, = 2.4magarcsec“^, and 

- /Ue,/ = 0.5 magarcsec^^. The resulting surface-brightness 
profiles and the measurements are shown in Fig. E.3a. 

The measurements agree well with the model that was 
convolved with PSFyom, throughout the full radial range. The 
scattered-light halo radius rno - 31" (rno - 41") with PSFyom 
(PSFkvi). M13a measure PSFmi 3 for the same telescope that is 
even fainter than PSFkti, but this PSF only seems to be correctly 
measured out to r ^ 30" and the times of the measurements 
that were used to derive PSFmi 3 are unknown (see Fig. A.l and 
Sect. 5.2.2). The measurements extend out to r = 70", which 
requires a PSF that extends to r = 77" (Sect. B). As with Mrk 
297, my study suggests that the measurements are scattered light, 
more accurate measurements of the PSF would reveal to which 
extent. Notably, the authors make an effort to carefully subtract 
the sky, which is why it should be possible to exclude this element 
as an explanation. 


1 achieve the same conclusion with the surface-brightness pro¬ 
files of BCG VCC 0001 that are presented by Meyer et al. (2014). 
These authors have chosen to only present the measurements of 
the host and not of the (possibly) bright starburst region, which 
makes an analysis of the origin of the excess light impossible 
(profiles are shown in the online-only Fig. E.5). 

The B - I colour profiles in Fig. E.3b show red ex¬ 
cess for radii r ^ 18". The red excess at r = 30" is 
A(B - 1) ^ 0.30magarcsec^^ with PSFyom and PSF,_om and 
about 0.24 mag arcsec^^ with PSFy 3 n, and PSF/^m. The same 
values at r = 40" are 0.79 mag arcsec^^ and 0.63 mag arcsec^^, 
respectively. 
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Fig. E.4. Surface-brightness profiles versus ra¬ 
dius for the alternative models of ESO 350- 
IG038 (Haro 11), compare with Fig. E.l. Pan¬ 
els, lines, and coloured bullets are as described 
in Fig. 2. The /-band profiles are offset by 
-1.5 magarcsec^^ from the B-band profiles for 
increased visibility. Grey and black bullets and 
error bars mark measured S-band (BO02) and 
V-band (MIO) values, respectively. The grey dot¬ 
ted (dashed) line shows the best-fit Sersic-type 
model (exponential fit to the disc in the range 
20<r<41")ofB002. 
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Fig. E.5. Surface-brightness profiles versus ra¬ 
dius for the BCG VCC 0001. Panels, lines, and 
coloured bullets are as described in Fig. 2. The 
i-band profiles are offset by -1.5 magarcsec^^ 
from the r-band profiles for increased visibil¬ 
ity. Black diamonds (grey; orange) diamonds, 
and error bars mark measured r-band (g-band; 
f-band) values of Meyer et al. (2014). 











